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Abstract. We consider a beam of two-level randomly excited atoms that 
pass one-by-one through a one-mode cavity. We show that in the case of 
an ideal cavity, i.e. no leaking of photons from the cavity, the pumping 
by the beam leads to an unlimited increase in the photon number in 
the cavity. We derive an expression for the mean photon number for 
all times. Taking into account leaking of the cavity, we prove that the 
mean photon number in the cavity stabilizes in time. The limiting state 
of the cavity in this case exists and it is independent of the initial state. 
We calculate the characteristic functional of this non-quasi-free non- 
equilibrium state. We also calculate the total energy variation in both 
the ideal and the open cavities as well as the entropy production in the 
ideal cavity. 
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1. Model and Results 
1.1. The Model 

Our model consists of a beam A of two-level randomly excited atoms that 
pass one-by-one through a photon cavity C. During the passage time r the 
atom in the cavity interacts with the cavity field. 

The cavity is a one-mode resonator described by a quantum harmonic 
oscillator with Hamiltonian Hq — eb*b <E) 11 acting on the Hilbert space He- 
Here b* and b stand for boson creation and annihilation operators with canon- 
ical commutation relations (CCR): [b , b*] = 11, [b,b] = [b*,b*] = 0. 

The beam of two-level atoms with energy levels and E > 0, is described 
by a chain Ha = ^2k>i HA k °f individual atoms Ak with Hamiltonian H^ k = 
1 £g> Erj k in the Hilbert space Ha = ®k>iH.A k - Here for any k > 1, H.A k = 
C 2 and the individual atomic operator r\ k ■= [<J Z + -0/2, where a z is the 
third Pauli matrix. The eigenvectors ipk : Vk^t = ^k an< ^ Vki^k = 0; are 
interpreted as the excited and the ground states of the atom, respectively. 

The initial state of the system is the product state of the cavity and the 
states of each individual atom: 



Here pc is the initial state of the cavity, which we assume to be normal, i.e., 
given by a density matrix pc G £i("Hc), the space of the trace-class operators 
on %c, and Pa '■= ®k>i Pk 1S the state of the beam. 

In this paper, we suppose that the atomic states {pk}k>i on the algebras 
{£i(%.4 fc )}fc>i are diagonal and identical (homogeneous beam), hence of the 
form 



The parameter < p < 1 denotes probability that atom in beam is in its 
excited state ip~£ . 

For simplicity we consider our model in the regime when at any moment 
only one atom is present in the cavity. Physically, this corresponds to a special 
tuning of the cavity size / and the interatomic distance d = I. Then the cavity- 
atom interaction has piecewise constant time-dependence, and we take it of 
the form: 



Ps ■= Pc®(£)pk ■ 



(1.1) 



k>l 




(1.2) 



K k (t) = Xp-i)r,feT)(*) A (6* +b)®n k . 



(1.3) 
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Here Xa{%) is characteristic function of the set A. 

Since the non-excited atom is not "visible" by the cavity (Wk(t)(u ® 
ip^) = for any u G dom(i?e)), the detuned case d > I, when there is 
at most one atom inside the cavity at the same time, is also described by 
a piecewise constant interaction. This situation can be handled by a small 
modification of our arguments, but we will not consider this situation further 
in this paper. 

Remark 1.1. It is convenient, albeit not compulsory, to stick with a quan- 
tum description of the atomic beam. If one restricts oneself to the atomic 
observables rjk, which generate a commutative algebra, one can also consider 



Ha and interaction (1.3 1 as matrix representation of the continuous-time 
Bernoulli process with time-unity r. Then the its piecewise constant random 
realisations 

V(t) =y2(i } k'X[(k-l)r,kr)( t ) T lk1pi)nA h > ( L4 ) 



fc>l 



are generated by random sequences of atomic operator rjk '■= {c z + -0/2 
eigenvectors {V't }fe>i with probabilities p and 1 — p for eigenvalues 1 and 



respectively. Here (• , -)u^ k 1S the scalar product in the space 'Hj\. k = C • . 

Returning back to the quantum electrodynamic origin of the interaction 
(1.31 one observes that it is completely elastic, since the atomic system does 
not evolve. The atom remains in the same state throughout its interaction 
with the photon field. This may be interpreted as the limit of rigid atoms 
that "kick" the cavity mode, see |FJMa| . 

The Hamiltonian for the entire system S acts on the space Hs '■= He <8> 
Ha-i an d is given by the sum of the Hamiltonians of the cavity and the atoms, 
and the interaction between them: 



H(t) = H c + J2(H Ak +K k (t)) (1-5) 

fe>i 

= eb*b(g)t + ^2t(g> Erik + Y, X[(fe-i)r,fcr) (f)(A (b* + b) ® %) . 



k>l k>l 



Notice that for the time t G [{n— 1)t, tit), n > 1, only the n-th atom interacts 
with the cavity and the system is autonomous. 



Projected on the time invariant sequences of {ip k }k>i the system (1.5) 



is a one-mode cavity Hamiltonian with random interaction driven by the 1-0 



Bernoulli process (1.4 1 



1.2. Hamiltonian Dynamics of the Ideal Cavity 



Let t € [{n— 1)t, tit). Then the Hamiltonian (1.5) for the n-th atom in the 
cavity takes the form H{t) = H n , where 



H n := eb*b® t + V t®En k + \{b* + b) (g) r?„ . (1.6) 



fc>i 



Although the atomic beam is infinite, since we will assume that the 
initial state is normal and since up to any finite time t only a finite number 
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of atoms have interacted with the cavity, we can describe the evolved system 
by normal states uis(') '■= Tr( • ps), which are defined by density matrices ps 
from the space of the trace-class operators €i(T-Lc <X> Ha)- Then the partial 
traces over Ha and over He : 

wc(-) := u s {- ® 1) = Tr c (- Tr APs ) , (1.7) 
u A {-) :=w s (l®.)=^(- Tr cPs) , 
define respectively the cavity and the beam states. 

Since below we mostly deal with normal states, we will use the terms 
'state' and 'density matrix' interchangeably, if this does not cause any con- 
fusion. 

We suppose that initially our system is in a product-state ps (t) | t _ Q = 
Pc ® PA °f the sub-systems C and A: 

£4(-)=Tr(- pc® pa) • (1-8) 

For any states pc on 2t("Hc) and p A on 2t("H^) the Hamiltonian dy- 
namics of the system is defined by (1.5|, or by the quantum time-dependent 
Liouvillian generator: 



L(t)(p s (t)) :=-i[H(t),p s (t)} 



(1.9) 



Then the state ps{t) '■— (pc ® PA)(t) of the total system at the time f is a 
solution of the non-autonomous Cauchy problem corresponding to Liouville 
differential equation 

j t p s (t) = L(t)( Ps (t)) , ps(t)\ t=0 = P c®PA- (1.10) 



We denote by ujg(-) := Tr( • ps(t)) the system time evolution due to (1.10) 
for the initial product state tos(') (1-8). 

Notice that in general (see e.g. |BrRolj ) the Hamiltonian evolution 
( 1.10 1 with time-dependent generator is a family of unitary mappings {Ut, s }s<t 
(called also evolution operators or propagators) 

. d 



i-U t , s = H{t)Uus , s<t, 



for s < r < t . 



with the composition rule: 

Ut,s = Ut, r U rjS , 

Here s,r,f6l and J7 tjt = 1. Then solution of ( 1.10 1 is 
Ps(t) - U t , sPs (s) U t -j =: T tia (p s (a)) 



(1.11) 
(1.12) 
(1.13) 



where by fll.l2| ) the mapping T M (.) = r t , r (T P ,,(-)). 

For our model with the tuned repeated interactions the form of the 
evolution operator (1.12) and the solution (1.13) are considerably simplified. 
Indeed, our system (1.5) is autonomous for each interval [(fc— l)r, fcr). Then 

)uvillian generator 

for t E [(& - l)r, fcr) , k > 1 , (1.14) 



by virtue of |L6| and (1.9) the Liouvillian generator 
L(t)(.) = L fc (.) = -i[H fc) 
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is piecewise constant (time-independent). For any t > one has the repre- 
sentation 



t := n{t)r + u(t) , n(t) = [t/r] and v(t) £ [0,r) , 



(1.15) 



where [x] denotes the integer part of x G R. Then by ( 1.13| ) and (1.14) the 
solution of the Cauchy problem (1.10 1 for t £ [(n — l)r, nr) takes the form: 



Ps{t)=T tfi {p c ® PA ) :=e v ^ 



By (1.121 and (1.14) the mapping T tj o is the composition 



k=l 



k=n—l 



of the one-step evolution maps defined as 



= e rLk 



and r t) ( n _i) T 



(1.16) 



(1.17) 



(1.18) 



Consequently, the evolution of the initial state of the system (1.8) can be 
expressed as 

u t s (-)= r b(-T t , (pc® PA)) • (1.19) 
The mathematical study of this kind of dynamics, for different types of 



repeated interaction K n {i) (1.5), was initiated in papers by Attal, Bruneau, 
Joye, Merkli, Pautrat, Pillet: |BJMfl Kfal KTT1 \KM IBJM21 iBPll IBJM3j . 
These works provide a rigorous framework for such physical phenomenon as 
the "one-atom maser", see e.g. [MWM, FJMbi. The important mathematical 
aspect of repeated interactions is a piecewise constant (random) Hamilton- 
ian dynamics like (1.16), which in certain limit of r — > may produce an 



effective quantum Markovian dynamics, which drives the total system to an 
asymptotic state. In our case, we will be able to obtain the exact asymptotics 
of the dynamics of our simple model directly, i.e., without taking a limit. 

In the next part of the present paper (Section [2]) we exploit specific 
structure of the Hamiltonian dynamics (1.16) and a special form of interaction 



( 1.3 ) to work out an effective Hamiltonian evolution of the perfectly isolated 



cavity C. In this case, the pumping of the cavity by the atomic beam leads to 
an unlimited growth of the number of photons. This means that the limiting 
state of the cavity will not be described by a density matrix. For this case, 
our results concern evolution of the photon-number expectation N(t) in the 



reduced by the partial trace (1.7), (1.19) time-dependent cavity state 



WcO) : = ® 1) = TrcG Pc(t)) for Pc {t) := Tr A (T tfi (p c ® p A )) . (1.20) 



Note that in our model, see (1.5) and (1.6), the atomic states pu do not 
evolve: 



[l®p k ,H(t)] = , k> 1 



(1.21) 



The form of the initial state ( 1.1 ) together with ( 1.16 ) and ( 1.20 ) deter- 



mine a discrete time evolution for the cavity state: := pc(t = nr), with 
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the recursive formula 

P ( c n) ~Tr APs (t) = Tr A [e TL " 



e™(pc®(X)p fc )] 



(1.22) 



k=l 



n-1 



=Tr^[e rL "{Tr^_ 1 ...T Ml 



■ e TL2 e TLl (p c ®Q<)Pk)}®Pn} 
fe=i 



Here we denoted the partial trace over the fc-th atom space H Ak by Tr^ fc (•). 

For any density matrix p £ £i("Hc) corresponding to a normal state on 
the operator algebra 2l(%c) we define the mapping C : p >-» £(p), by 



£(p) := Tr^ fc (e ri * (p ® p fc )) = Ti Ufc [e 



-irHu 



(p® Pk)e 



irH k 



(1.23) 



Here the last equality is due to (1.14 1. Note that the mapping ( 1.23[ ) does 
not depend on k > 1, since the atomic states {pfc}fc>i are homogeneous ( 1.2 1. 
Then the cavity state at t = kr is defined by the fc-th power of C: 



p^=c( P r i, )=^(p C ) 



(1.24) 



Therefore, by (1.161, (1.221 and (1.24) one obtains that for any time t = 
(n — l)r + v(t), where v{t) € [0, r), the cavity state is 

Pc(t) = ^„[e"« L »(£"- 1 (p c ) ® p„)] . (1.25) 

Our first result concerns the evolution of the expectation value N(t) of 
the photon-number operator N := b*b in the cavity at the time t: 

N(t) := 4(6*6) - Tr c (6*6 Pc {t)) . (1.26) 



For t = tit the state of the cavity can be expressed using (1.24). Then (1.26) 
yields 

N(nr) = Tv c (b*b C n {p c )) . (1.27) 
In the theorem below we suppose that the initial cavity state w^| t= o(-) = <jJc{') 
is gauge invariant, i.e. e laN pce~ %aN = pc- 

Theorem 1.2. Let pc be a gauge-invariant state. Then for a homogeneous 
beam the expectation value ( 1.27) of the photon-number operator in the cavity 
at the moment t — nr is equal to 



2A 



2A 2 



N(t) = N(0) + np(l - p) -^(1- cos er) +p 2 ^-(1- cos ner) . (1.28) 

e z e z 

If for the initial state pc one takes in the theorem the Gibbs state for 
photons at the inverse temperature /3: 



4 



-finb'b 



/Tr c e 



-Peb'b 



(1.29) 



then the the first term in (1.28) is simply JV(0) = {e^ - 
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If the initial cavity state pc,r,4> is not gauge-invariant with the breaking 



gauge-invariance parameter: r e 



of (1.28) one gets by Lemma 2.3 
N(t) 



w c,r,<ji(&) — Tr^(6 PCr,s)i then instead 



where ux 



2\r . 
p [cos ( 



cos(ner 



(1.30) 



■C,r,*( N ) U=0 



N(0) > r 2 . The same lemma yields also the positivity 



of (1.30) for any time t > 0. 

Remark 1.3. If all atoms in the beam are in the ground-state 



)fc>iV>fc , 



one 



has: (Er)k)ip k = 0, and pk | P =o = (V^ r^A.^Pk ' s tne projection operator. 



Then by (|l.28| and (|l.30| w^ r ^(JV) 



Since the mean cavity 



energy is e(t) = Trc {He pc (t)) = e -^(i). Therefore, there is no energy transfer 
from the beam to the cavity when p — 0. Whereas in the other extreme case, 
when all atoms are excited Pk{f]k) = L we obtain 

2A 2 



UJ, 



C,r, 



=0J. 



C,r,<f>\ 

2\r , 



=o(iV) 

t - 

; — cos(ner 



1 



cos ner 



(1.31) 



*)] 



In case of the gauge-invariant initial state (r = 0) the mean value of the 
photon number (1.31) satisfies the estimate N(t) > N(0) and it is oscillat- 
ing between initial value N(0) and N(0) + 4A 2 /e 2 with the cavity resonant 
frequency e. 

Remark 1.4. These Rabi oscillations (see e.g. |F JMb| ) are a simple conse- 
quence of non-trivial Heisenberg time evolution of the cavity boson operators 
generated by Hamiltonian (1.6). For example, if t € [0, t), i.e. n = 1, then we 
obtain (cf. Lemma 4.1): 



t : bi 



itH_ 



= e tL "=i(6® 11) = e~ iet b® 1 
Here evolution maps on the algebra 2l(%c ) 



-UHi 



n 



T, 



t,(n— l)r 



'(ft® l)e 
- 11® -?7i(l 

'St(^): 
(n-l)r- 



i/(t) , 



(1.32) 



(1.33) 



are adjoint to (1.17), (1.18) by duality with respect to the state (1.8), see 
Appendix A. 2. Here L* k denotes operator, which is adjoint to the Liouvillian 
generator ( fLli| . Then (Tig) and ( |LTf| , pT33] ) yield 

n-l 



c4(-)=Tr(T t y.) pe®P,0 



(1.34) 



fe=i 



Remark 1.5. The Hamiltonian evolution (1.32) breaks down the gauge invari- 
ance of the initial state wc(-) = Trc( • pc)- Indeed, by (1.7) and (1.8) one gets 
for the initial gauge-invariant cavity state: wc(b) 



|L20| and ( |1.34[ ) we have 
4(6) 



ws(&(g> 11) = 0. Since by 
w|(6®l)=Tr(T* (6®l) p c ® , (1-35) 
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(L2j) and (|1.32[) imply that for t G [0, r) 



wc(&) =J»-(e" 



1) 



(1.36) 



This property of dynamics is due to interaction (1.6 1, and it has also a non- 
trivial impact for the open cavity evolution. 

1.3. Quantum Dynamics of Open Cavity 

To make a contact of our model with a more physically realistic situation, we 
consider an open cavity. This allows the photons to leak out of the cavity, but 
also to diffuse in from the environment. In Theorems 11.61 and 11.81 we show 
that if the leaking rate is greater than the rate the environmental pumping, 
the photon number in the cavity stabilizes at a finite mean value. 

We consider this case in the framework of Kossakowski-Lindblad exten- 
sion of the Hamiltonian Dynamics to irreversible Quantum Dynamics (see 
Appendix A.l and [AlFal iAjPIIj . ) with time-dependent generator: 



L a (t)(p s (t)) :=-i[H(t),ps(t)}+a- &®lp s (i) b* ® 11 - 



-{bb* <g> l,p s (t)} 



l Ps (t) b®l-^{b*b®l )Ps (t)} , 



(1.37) 



for the complete positive evolution of total system, when the initial state 
is ps{t — 0) — pc ® PA (1-8 1. For a T > the cr-part of this generator 



(cf (1.9)) corresponds to the non-Hamiltonian part of dynamics. Here cr_ 



describes the rate of the photons leaking out of the open cavity into the 
environment, whereas o~ + corresponds to the cavity pumping rate due to the 
photon infiltration from the environment. It has to be distinguished from the 
pumping mechanism due to the interaction with the atomic beam. 



Note that similar to (1.10) the evolution of the state is defined by the 



solution of the non-autonomous Cauchy problem corresponding to the time- 



dependent generator (1.37). The proof of existence of this solution is a non- 



trivial problem (Appendix A.l), see for example [NVZj for the case of lattice 
systems and bounded generators. For unbounded generators and for a general 



setting the proof that the Kossakowski-Lindblad generator in the form ( 1.37) 
corresponds to a properly defined continuous evolution is more involved NZ, 
VWZ . A separate problem is to prove that this map is trace-preserving and 
verifies the property of complete positivity, which are indispensable for correct 
description of the open system evolution, see Appendix A.l and the references 
there. 

To avoid these complications we consider the case of the tuned repeated 
interaction, when the Hamiltonian dynamics is piecewise autonomous for each 



interval [(fc — l)r, fcr). Then for t € [(k — l)r, kr) the generator ( 1.37 1 has the 
form: 



L CT , fe (0 :=-i[H k ,-]+a- b ® 1 (•) b* ® 1 - 
+ 0-+ b* ® 11 (•) b® 1- ^-{b*b( 



2 



{bb* 



5 A,-} 
k > 1 



(1.38) 
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a.k 



— Lk, see (1.14). Hence, the solution of the non-autonomous 



and L, 

Cauchy problem 

3lAS,<r(*) = La{t)(p S ,<r(t)) , PS,a(t)\ t=Q = PC <8> PA 



(1.39) 



for the piecewise constant generators (1.38), has the same form as for the 
ideal cavity, a T = ( |1.16| , fll.l7| ): 

p s At)=TZ {pc®PA)-= (1-40) 



Here t = (n — 1)t + f (t) and the mapping X^q is the composition of the 
one-step non-Hamiltonian evolution maps defined by (1.38): 

tL„ 



rp(J rp(T 

1 k ■— 1 kT : (k-l)i 



and T° 



= e v{t)L a 



(1.41) 



-t,(n-l)T 

By duality with respect to the initial state ujg(-) = Tr(- pc S3 p.4) one 
can now define the adjoint evolution mapping {(T° )*} t >o by the relation 

^JA) = Tr(A T° (p c ® p^)) = wg^o) V)) , (1.42) 

for any A € 2l("Hc ® %.a)- Then similar to the nonleaky case ( |1.34 1 we obtain 
for any t = (n — l)r + v(t) that 



n-l 



<4,a0 - Tr ((T° o y(-) Pc ® AA) , TO* = [] e 



(1.43) 



fc=l 



Here {L* ak }k>i are generators, which are adjoint to ( |1.38 ). 

Since the restriction of (1.40) to the dynamics of a cavity state p € 
£i(Hc) is t ne partial trace over beam states, the corresponding discrete evo- 
lution mappings (1.23) have to be modified for the open cavity as follows: 

C*{p) :=Tr Ak (e TL ^(p®p k )) • (1-44) 

The mapping (1.44) does not depend on k > 1, since the atomic states 
{pk}k>i are homogeneous (1.2). If pc := Pc.a(t) |t=o is initial state of the 
open cavity, then similar to (1.25) we obtain by (1.44) for pc.ait) at the 
moment t — (n — l)r + v(t): 

Pc At) = Tr An [e v ^AC~ 1 (Pe)^Pn)} ■ (1-45) 
By ( 1.45 1 and (1.44 ) we obtain for the time-dependent open cavity state 
: = <4,<r( ■ ® 1) = Tr c (- p c , ff (i)) , (1.46) 
where u)* s „(■) := Tr( • T t ff (p c ® p^)) by Ool. We also define 



w o(') : = Tr c( ' Pc) and w C)CT (-) := lim u e (•) 



(1.47) 

To study the infinite-time limit wc. (7 (-), we consider the functional 

we lt r(W(0) = lim ^cAW(O) , (1.48) 
generated by the Weyl operators on He'- 

W(() = e^ b+Cb " ) , C € C . (1.49) 
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Notice that convergence ( 1.48 ) on the family of the Weyl operators guarantees 
the weak limit [BrRolj of the states when t — > oo, see Appendices A. 2 

and A. 3. The following theorem is our first result about the open cavity. 

Theorem 1.6. Let cr+ > and cr_ — a + > 0. Then for any gauge- invariant 
initial cavity state pc and for a homogenous atomic beam with parameter 



p = Tru Ari (r) n Pn), the limiting cavity state 

"Cat) : = , lim w C,<x(') 



(1.50) 



exists and it does not dependent on pc- Here the limit means trace-norm 



convergence of the sequence (1.4-5) to a density matrix pc.o 

lim \\p c ,a - Pc,a(t)\\i = 0, 



(1.51) 



functional (1-48) is 



where the norm \\ ■ ||i denotes the trace-norm on the space of trace-class 
operators Ci('Hc); see (5.31), Appendix A. 2. The explicit form of the limiting 

(1.52) 



x E{ pexp ^G (1 

where p := ie + (er_ — cr + )/2. 



-flT\ 



This result is obviously different from the case: er_ = cr + = 0, when 
there is no regular limiting state, because the number of photons in the ideal 



cavity cavity unboundedly increases with time, see Theorem 1.2 



Corollary 1.7. Theorem 1.6 implies that wc.o-O) is a regular, normal (see 



Appendix A. 5) and, in general, non-gauge-invariant state. One also sees that 
it is not quasi-free for < p < 1, but it obviously does for p = Q or p = 1. 

To verify these properties notice that by the Araki-Segal theorem (see 



Theorem 5.2 1, any regular state over CCR('Hc) is uniquely defined by the 
characteristic functional H> ujc^(W(())}^ e(: .. To check the conditions of 



this theorem we express the infinite product in (1.52) as a uniformly con- 



verging (by estimate (3.27)) infinite sum of terms that have a generic form: 



'+ P n (l-p) m e 



where 



r fc (C) := 2 Im ([A(l - e^eT^C]/^' 



(1.53) 



(1.54) 



By virtue of (5.55), apart from the normalization, expression (1.53) is noth- 



ing but the characteristic function of a non- gauge-invariant (see (1.54)) and 
quasi-free state on CCR("Hc), which trivially verifies the Araki-Segal theo- 



rem. As a consequence, (1.52) defines a normal state since it is a convergent 



infinite convex combination of normal quasi-free states but, when < p < 1, 



the infinite combination of quasi-free states in (1.52) is not quasi-free. 
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Let N a (t) be the expectation value of the photon- number operator b*b 
in the open cavity at the time t as: 

N a (t) := Lu c Jb*b) = Tr c (b*b p c , CT (i)) . (1.55) 

For the open cavity the result corresponding to the asymptotic be- 
haviour of the photon number in the cavity takes the form. 

Theorem 1.8. Let cr_ — <r + > 0. Then for an arbitrary initial gauge-invariant 
cavity state pc such that the initial mean-value of the photon number in the 
cavity is bounded : 



u C,a 



(b*b) 



Tv c {b*b p c ) < oo 



we obtain that the expected number of photos att = (n— l)r 
form: 

N a (t) = e-^-^N^O) 



(1.56) 
v(t) has the 

(1.57) 



A 2 

p ri2 ( 



- ff +) T (l-e" T )(l 



1 



-(<r_-er + )t 



,2A 2 1 



-(er_-<r + )t 



2 2A 2 n 

p w {1 - e 



:(1 



-(<>•- 



g-(CT_-CT + ) 



' 1 — e -(<T_-0- + )T 

ff +5 T / 2 coser) 



"(1 



(<T_-<T + )t\ 



TJere /i := (cr_ — o~+)/2 + ie. TVie /irmi of the expected number of photons in 
the cavity is 



ucAb*b) := lim uj c a (b*b) 

t— too ' 



(1.58) 



= p(l-p) 



2A 2 1 - e -( CT -- a +) T / 2 coseT 



1 — e -(C--f+)T 

the limit 

p(l -j>) 



+ p(2p-l) 1 



<T_ — (J 4 



Remark 1.9. Note that the limit (1.58) satisfies the estimate from above: 
2A 2 



M 2 i 



-(<t_-ct + )t 



p(2 P - 1; 



A^ 

\iA 



as well as the estimate from below: 
w c , CT (6*6)> ' 



(1.59) 



(1.60) 



| M |2 I + e (a_-<r + )r ff _ _ CT+ ' 

There are several instructive cases that one can observe as corollaries of the 
above theorem: 

• Let p — (or A = 0), i.e. no pumping by the atomic beam. Then the 
limiting value ( 1.58 1, ( 1.60 1 is <7 + /(<7_ — o+), which coincides with the 
formula (5.42| in Appendix A. 3. If in this case one put rr_ > and 
consider 0+ — > 0, then the limit (1.58) is equal to zero. This simply 
means that in the absence of atomic pumping the leaky cavity relaxes 
to the vacuum state empty of photons, or formally to the Gibbs state 
(5.44) with zero temperature /3 cav = +00. 
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When A ^ and p = 1, still there is no (unlimited) pumping by the 
atomic beam but only the photon Rabi oscillations ( 1.28 ). Then for the 
open cavity er_ > c+ > the limit ( 1.58| ) and (1.59) give 



A 2 



0"+ 



(1.61) 



for any leaking cr_ > 0, including the limit of the ideal cavity when 
cr_ — > 0. Note that in the ideal cavity for the case p = 1 the mean-value 



of photons (1.28) is bounded, but oscillating. In other words, the limits: 
t — > oo and ct_ — > do not commute. 



If < p < 1 and er + = 0, then for the limit of ideal cavity ( 1.58 ) yields 

2A 2 ,. (1-coser) 



lim u>c,a{b*b) = p(l - p) 

7-— >0 



lim 



>o 1 



(1.62) 



Hence, for the non-resonant case er ^ 2its, where s £ Z, this limit is 
infinite, i.e. corresponding to the conclusion of the Theorem 1.2 about 



unlimited pumping of the ideal cavity. Indeed, ( 1.57) for <r + = implies 
lim 4=" T (6*6)| ff+=0 - N(nr) . (1.63) 

a — — ^0 



Here N(nr) coincides with ( |1 .28 1, which diverges for the non-resonant 
case as np(l — p) 2A 2 (1 — cos er)/e 2 |„_ i . 00 , cf. (1.62). 
When a + > and <j_ — > a + , the limit ( |1 .58 1 yields 



lim ujc a(b*b) = +oo 

--HT+ 



(1.64) 



It means that if the leaking and the environmental pumping have the 
same rate, the limiting state corresponds to the infinite temperature 



state, see (1.52) and (5.43), (5.44). On the Weyl operators this state is 



given by the Kronecker delta-functional 



i if c = o, 

if C 7^0. 



Since the Kronecker characteristic functional is not continuous, the cor- 
responding state is not regular. Consequently, the Araki-Segal Theorem 
|5.2| is not applicable. 
Note that for er + > and tr_ 
number of photos at t = nr 



<7+ by (1.57) one gets for the expected 



lim 



«6*&) 



iV(nr) + n t<j + 



(1.65) 



Hence, in this case the pumping by the random non-resonant atomic 



beam (1.28) and by the environmental pumping due to ct + > 0, give the 



same linear rate for increasing of the mean number of photons in the 
cavity. Consequently, for t = nr — > oo the infinite photon number cavity 
state coincides with the infinite-temperature state that we discussed 



above, cf. (5.44) 
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This concludes the description of our main results. The rest of the paper 
is organized as follows. In Section[2]we use the specific properties of our model 
to diagonalise it, which is the key for the further analysis. Then we give the 
proof of Theorem \1.2\ for the Hamiltonian dynamics of the ideal cavity. 

In Section [3] we present our results for the case of the leaking cavity 
described by the Kossakowski-Lindblad irreversible quantum dynamics, i.e., 
Theorems [T6] and [L8l 

The results concerning Energy-Entropy relations are presented in Sec- 
tion |2J There we calculate the energy variation (Theorem |4.2| and Theorem 
|4.8| and obtain a formula for the entropy production (formula ( |4.47 1) for the 
ideal cavity. 

The Section [5] is reserved for comments, remarks and open problems. 
For the reader convenience we collect in an Appendix certain results 
and definitions necessary for the main text. 



2. Hamiltonian Dynamics: The Ideal Cavity 

Here we consider the case of the ideal cavity, i.e. cr_ = <r+ = 0. In this 
case, the discrete evolution map £ CT= o — C is given by (1.23). Recall that 
p = Tr(rj k p k ) for the homogeneous atomic states {p k } k >i in the beam. 

Our first result co ncern s the expectation of the photon- number operator 
N = b*b in the cavity (1.26). For t = nr this expectation involves the calcu- 
lation of C n (p) ( 1.27). Instead, we use the n-th power of the adjoint operator 
C* defined by duality with respect to the cavity state ojc(-) — Trc( ■ pc), see 



(1.271 and Remark 2.2 



Let Sk be ^-isomorphism (unitary shift) on the algebra 2t(Hc) ( X>2l(H^i fc ) 
defined by 



S fe (-):=e^C)e" tVfc , V k := \(b* - b) ® Vk /ie 



,-tV. 



Since rjfc = rjk and ® %) = 1 ® whereas 
S k {b®l) = b®l-t® X r lk , S k (b*®l) = , 



1 



„ A 

1® -r] k 
e 



the transformation (2.1) of the Hamiltonian |L6j gives 

A 2 



H k := S k (H k ) = eb*b®l+l®(E 



^2 ^®E-q s 

s>l:s^ k 



(2.1) 



(2.2) 



(2.3) 



Notice that the dynamics generated by H k (or by H k ) leaves the atomic 
operator r\ k = (a z + 1)/2 invariant 



(1® 7] k )e 



-irH k 



1® Vk 



Similarly to the unitary shift (2.1), we define on 2l("Hc) the ^-isomorphism: 



S(-) 



(•)e- 



V := \{b* - b)/ie 



(2.4) 
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Lemma 2.1. For any state p on 2l(%c) the one-step evolution mapping C 



1.23) has the form: 



C{ P ) = p S- 1 {e-' lTeb " b S{p)e lTeb ' b ) + (l-p) e - iTeb * b pe ireb * b . (2.5) 

Here p = Tr^ Afc (rj k Pk) is the probability to find the k-th atom in the excited 
state with energy E. 



Proof. Using the shift transformation (2.1 1 and the cyclicity of trace, one can 
re- write ( 1.23 ) as 

C{p) = Tv Ak [Sz\e-^S k (p®p k )e^)] . (2.6) 

Since the operator ra& = [a z + I)/2 is idempotent (r/^ = rjk) and since it 
commutes with pk (1.2), the combination of expansions of (2.1|, ( 2.4 1 with 
definitions of 14 and V gives 

Sk{p ® Pk) = S(p)® r\ k pk + P ® (I - m)pk , 

S^ip® Pk) = 5 -1 G°) ® VkPk + P®{! - Vk)Pk ■ 



(2.7) 
(2.8) 



Therefore, plugging ( |2.7[ ) into (2.6) we obtain 

C{p) ^r Ak [S^{e-"^{S{p)® m p k )e^)] 

+ Tr Ak [S7\e- iT ^ (p ® (/ - Vk)Pk)e iTS *)) 



Now, the diagonal form (2.3) of the Hamiltonian H k and (2.8) for S k , imply: 
C(p)=Tr Ak [S- 1 (e- i ™»' b S(p)e i 



eb*b\ 



+ Tr Ak [e-^ b b S(p)e 
+ Tr Ak {S- 1 (e~ lTeb * b 



® VkPk\ 

" lb * b ®(I-r)k)VkPk\ 

treb*b\ 



pe 



)®r}k(I ~Vk)Pk] 



Tr Ak [e-^ b " b pe^ b ®(I-n k )pk] 



Since (/ - r] k )rik = and p = Tr Ak [1 ® n k pk], one gets (2.5). 



□ 



Remark 2.2. (a) To calcu late t he expectation of the photon-number operator 
N = b*b at t = tit using ( 1.27), we would need to find the action of the ra-th 
power C n {p) of the operator (2.6). 

We show that in fact it is easier to calculate this mean value using the ra- 
th power of the adjoint (or dual) mapping C* . For any bounded operator 
A £ B{T-Lc) and p £ £i(Hc), it is defined by relation 

Tr e (C*(A)p):=Tre(A£(p)) , (2.9) 

see Appendix A. 2 for discussion and details. 

(b) Using invariance of the atomics states (1.21) one can obtain explicit ex- 
pression for the one-step dual mapping C* . Indeed, by (1.21) and by ( |1.23[ ) 

Tv c {A C{p)) =TT Hc ^ HAk {{A®t) e-" H *{p®t){l® Pk y TH *} (2.10) 

=^H c ®H Ak {(1 ® Pk)e- iTH " (A ® ty THk (p ® 1)} 



=Tr 



{e lTH «{A® p k ) e- lTH *{p® 1)} 
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with (2.10) yield expression for the one-step mapping 

£*(A) = Tr Ak {e tTHk {A®p k )e- lTHk } 



where we used cyclicity of the full trace Tr-^®?^ . Hence, (2.9) together 

g 

(2.11) 



which according to (1.23) is independent of k > 1 



(c) Let p £ Ci('He) be density matrix such that Trc("P(6, b*) C(p)) < oo for 
any polynomial V(b,b*). Then one can extend definition ( |2.9[ ) of C* to the 
class of unbounded observables, which contains V(b, b*) and in particular the 
photon-number operator b*b. The advantage to use the adjoint mapping C* 
is that its consecutive applications do not increase the degree of polynomials 
in variables b and b* . 



Now, applying to (2.11) the same line of reasoning as in the proof of 



Lemma [2.1 1 we find explicit expression for the adjoint one-step mapping 



C*{A) = V S^O* b S(A) 



eb'b 



)+ (I-P)e 1 



eb'b 



Ae 



tb'b 



(2.12) 



Note that alternatively one can obtain (2.12) using (2.5) and definition fl2.9|. 



Lemma 2.3. For A = b*b and for the adjoint operator C* defined by (2.12) 
we obtain: 



(£*) n (b*b) = b*b + p ^[(l-e meT ) b* + (l-e 

2A 2 2 
+ np(l - p) -^-(1 - cos ex) + p 2 - 



i£T ) b] (2.13) 
-(1 — cos Tier) . 



Proof. Since ^-isomorphism ( |2.4[ ) is a shift transformation, one gets 

S(b*b) = (b* - X/e)(b- X/e) . 
The CCR relations for b* and b yield: 

e r T ,b'b h * e -rr,b*b = jT % * ^ ^reb'b^-ireb'b = e ~ir % 

Consequently, 

S - l {e lTtb '* b S{b*b)e- lTeb ' b ) = fiT 1 ((e ,r6 &* - X/e)(e' lTe b - X/e)) (2.14) 
= (e tTe b* - (1 - e lTt )X/e){e- lTt b - (1 - e~ tTe )X/e). 



Hence, (2.14) implies for (2.12) 



C*{b*b) = p(e leT b* - (1 - e leT )X/e){e~ leT b - (1 - er ltT )X/e) 
+ (l-p) b*b 



(2.15) 



A A 2A 2 

b*b+p-(l- e leT )b* +p-(l- er l ")b + p^-(l - cos er) 



which coincides with expression (2.13) for n = 1. If we insert in (2.12) A = b* , 



and then A — b, we obtain correspondingly: 

C*{b*) =e ltT b* -p(l-e leT )X/e , 
C* (b) = e- ltT b - p(l - e" ieT )A/e . 



(2.16) 
(2.17) 
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Now one can use induction. Since (C*) n (b*b) = C*((C*) n ~ 1 {b*b)), we can 
apply ( |2.15| )-( |2.17| to \2.lty for n- 1 to check the formula \2.\i\ for the ra-th 
power. □ 



Remark 2.4. As we indicated in Remark |1.5| the evolution C breaks the gauge 
invariance of the initial cavity state pc ■ For n = 1 the gauge breaking param- 
eters r n= \ and <j) n =i ar e defined by (1.36), or by (2.9), (2.17): 



L>Z(b)=Tr e (C*(b)p c )= P ~(e 



1) 



r± e 



(2.18) 



Formulae (2.161,(2.17) allow to calculate r n and </>„ for any n > 1. Iterating 
them one obtains 



(£*)"(6) 



P (l 



e meT )A/e , 



b — p(l — e 



")A/e 



(2.19) 
(2.20) 



Therefore, the gauge breaking parameters r n and 4> n are defined by equation: 



L>r(b) = Trc((£*) n (b)pc)=P*(e 



-1) =:r n e**- 



(2.21) 



Proof, (of Theorem 1.2) To find the mean-value of the number of photons 
in the cavity at the time t = nr, we calculate the expectation of (2.13) in 
the initial cavity state ojc{-) = Trc(- pc), (2.9). Since pc is supposed to be 
gauge- invariant, this yields 



N{t) =Tr c (6*6 C n (p c )) = Tr c ((£* )"(6*6) Pc ) 

2A 2 2A 2 
=A^(0) + np(l — p) (1 — coser) +p 2 ^-(l — cosner) 



which coincides with (1.28) for t = nr. 



(2.22) 



□ 



Corollary 2.5. If the initial cavity state pc.r,4> is not gauge-invariant, then 
(|1.30|) and (|2.22[) give for i = nr 



2A 2 2A 2 
AT(t) =A^(0) +np(l-p) — (1 - cos er) + p 2 (1 - cos ner) 

2Xr r . 
+p [cos — cos(ner — <p)\ . 



(2.23) 



Remark 2.6. Theorem 1 1 . 2| implies that the pumping effect is non-trivial (i.e. 
N(t) is increasing to infinity) only for the beam of randomly exited atoms, 
i.e. < p < 1. There is no unlimited pumping of cavity by the non-exited, 
p = 0, or totaly exited, p = 1, atomic beams. Note that by (1.5) the mean 
value of the photon cavity energy is defined by £c(t) '■— eN(t). Hence, its 
unlimited increasing is due to the kicking by randomly excited chain of rigid 
atoms, which are pushed through the cavity. 
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3. Irreversible Quantum Dynamics: The Open Cavity 

In a real physical cavity it is always possible for photons to leak out of the 
cavity as well as to diffuse in from the environment. This realistic situation 
is modeled by what is known as an open cavity. In the present section we 
include these effects, via the Kossakowski-Lindblad extension (1.371 of the 
Hamiltonian dynamics ( 1.5 1 with generator ( 1.9 ). Below we consider the case 
of the open cavity with dominating photon leaking: cr_ > cr+ > 0, although 
our results allow also to study certain limiting cases, see Remark|l.9| 



Similar to the ideal cavity, see (2.5| in Lemma 2.1 we first find a new 
expression for the one-step evolution mapping in the case of the open cavity 
(L44I. 



Lemma 3.1. For any state p € £i(Hc) the one-step evolution mapping ( 1.44) 
for the open cavity has the form: 



C a {p) =pS- 1 {e TL ^{S{p))) + (1 -p) e T ^(p). 



Here S is defined by (2.4) and L\ M acts on 2t("Hc) as follows 

A/e) 



L\, a (p) ■= - ii^b.p] 

+ a-(b-X/e)p(b 

+ a + (b* -X/e)p(b-X/e) 

with L „ := La=o,ct- 



(3.1) 
(3.2) 



(T — 

~2~ 
°± 
2 



{{b-X/e)(b*-X/e),p} 
{{b*-X/e)(b-X/e),p}, 



Proof. Using ^-isomorphisms (|2.1|,(|2.4|) and equations (2. 7), (2. 8) we define 

(3.3) 



instead of generator (1.38) the following operator: 



La,k(p® Pk) ■= S k (L (Ttk (S k 1 (p® p k )) 
= -i[eb*b ® t + 1 (g) (E - X 2 /e)r, k ,p ® {n k p k + ( I - T] k )p k )} 

+ (a_(b-X/e)p(b* - A/e) - ^{(6- A/e) (6* - A/e),p}) ® mPk 
+ (a + (b* - X/e)p(b - A/e) - ^±{(b* - X/e)(b - A/e), p}) ® mPk 

a-bpb* - ^-{bb*,p}^ ®{I-r, k )p k + (a + b*pb- ^-{b*b,p}) <g> (I -n k )p k 



L\Ap) ® VkPk + L 0>a (p) <g) (/ - T] k )p k . 



With help of (3.3) and (2. 7), (2. 8) we obtain the representation 

e TL °- k (p <E> p k ) = S7 1 {e rt """ (S k (p <g> p k ))) 



= S k 1 (e^° 
= S k 1 (e^ 



k V c \^k\h"<y Hk) ) ) (3.4) 

(S(p) <g> T] k p k + p®(I - T] k )p k )) 

(S(p) <8 VkPk)) + S k \e TLa - (p) ® (I - r, k )p k ) 



S-^e^iSip))) <g> VkPk + e TL °«-{p) <g> (/- n k )p k 
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Therefore, the one-step mapping (1.44) takes the form 
C a (p) = Tr Ah (e rL -*(p®p k )) 
= p S-He rLx " (S(p))) + (1-P) e TLo - (p) 



which is claimed by Lemma in (3.1 1. 



□ 



Corollary 3.2. To define the adjoint mapping £* we use duality relation (2.9) 
for (3.1) and any bounded operator A G B(%c)' 



Tr c (CUA)p) :=Tt c (A C„{p)) 



(3.5) 



Then by definitions (2.4), (3.5) and by explicit form (3.1) of the operator 
C a (p) one gets that 



Cl{A)=pS-\{e^ 

If for adjoint operators (e TLx -<*) 
adjoint generators by 

e Ti U : = ( e riA .° 



■(S(A))) + (l-p), 



)*(A) 



(3.6) 



and (e T "'")* we define the corresponding 



and e T °-° := (e 



Tin 



then (3.1), (3. 2) and (3.5)-(3.7) allow to find them explicitly: 

L\ a {A) = « A] + ^(b*- X/e)[A, b] + U —[b\A]{b X/e) 

+ ^(b- X/e)[A,b*] + ^±[b,A}(b* - X/e) , 
L* 0t(T (A) = i[eb*b,A] + ^b*[A,b] + ^[b*,A]b 



(3.7) 
(3.8) 

(3.9) 



b[A, b* 



2 2 
by straightforward calculations. 



[b,A}b* 



Remark 3.3. As we indicated in Remark 2.2 (c), one can extend the ad- 
joint one-step mapping £* to the algebra generated by polynomials in the 
annihilation and creation operators. 



Remark 3.4. Duality allows us to define by the relation 
Tr(S* k (A) p® Pk ) :=Tr(AS k (p® Pk )) , 



(3.10) 



the operator S k with domain consisting of operators A — A <£> a E 2l(%e) ® 
St(H^ t t . ). By definition (2.1 1 an d by c yclicity of th e trac e, one obtains from 
( 3T0| that S% = S^ 1 , c£ ( 2~7| , |2~j| . Then Q and \3A0\ yield for the 
adjoint of the operator defined in (3.3) the following expression 

L* 



J a,k\ 



S k (Ll !k (S^(-))) 



(3.11) 



Here L* k is defined in (1.43). It can be explicitly calculated on operators 
A® a with help of ( |1.38[ )3mce ( pO} and ( |3.11| imply 



K, k (A ® o) = I-A j(T (i4) <8> %a + ® (J - »7fc)a 



(3.12) 
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this equation gives an alternative way to establish (3.8) and (3.9). Similar 
to (3.4), the representation (3.12) indicates that for each k > 1 the mapping 

{e TL " k } is reducible by two sub-algebras of $l(%c) ® 2l("H^i fc ): 



■■ k : A (g) rj k a \-> e rL ^- (A) ® r\ k a , 



n rL* 



" : A ® (I - ?7 fe )a h-> e Ti o,„ (A) ® (I - ?7 fe )a 



(3.13) 
(3.14) 



Proof, (of Theorem 1.6) Using the Baker-Campbell-Hausdorff formula one 
can rewrite the Weyl operator in the form 



W(C) = e l(cb+cb * )/v/ 2 = ^C&Vv^Cb/v^g-ICI 2 ^ 
From the CCR we find that for any f3,-f E C 

e" 6 *& = (b - /3)e' 9fc * and e~< b b* = {b* + 7 )e 76 . 



By (3.15) and (3.16) the action of L* x a is 



(3.15) 
(3.16) 



i A(er_ — a 



V2 



2c 



* ll {W{Qb--j=pCb*W(Q 

-(( + 0^(0 ■ 



^-\C\ 2 W(0 (3.17) 



Recall that fi — ie+ (er_ — cr+)/2 . Therefore, the mapping generated by ( 3.8 ) 
(and correspondingly by (|3.9[)): 



7a, 



(3.18) 



is quasi-free |DVV1[ IDVV2| . In particular, there exist functions ((t) and 
Qa,t(C) such that 



7A ) r(W(C)) = e- n ^^W(C(T)). 



(3.19) 



See Appendix A. 4 for definitions and details. 

To check this claim we calculate explicitly £(t) and fl\^ T (() using dif- 
ferential equation for (3.19): 



dr 



(3.20) 



By (3.17) the right-hand side of (3.20) is given by 



>( 



^,AC) ( _^_ K{T)w{C{T))b _ A K(t)6 ^ (((t)) 



2 lC(r)| a W(C(r)) + ^^ 



2c 



(C(r)+C(r))W(C(r))). (3.21) 
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The derivative of the r-dependent Weyl operator W(((t)) is calculated using 
Baker-Campbell- Hausdorff formula (3.151 

dW(((T)) 
dr 

Therefore, 7a,t(W / (C)) satisfies the following differential equation 



dr 

-QxAO ( _L^W 



V2 dr ^(CW) + ^^(C(r))6 



(3.22) 



Due to ( |3.20[ ), we can match the right-hand side of (3.22) with (3.21) 
and obtain the following system of differential equations for the functions 
C(r) and n A)T (C): 

dC(r) 
dr 

and 



-MC(r) 



gMQ l rf(C(r)C(r)) tx + a t A(cr_ ~ °"+) 

+ T lc(r)l £ — 



The solution to the first differential equation is 

C(r) = e-^C, 
and using it in the second equation we find 



(C(r) + C(r)). 

(3.23) 



d»A,r(C) 



,-(er_-<7+)T|^|2 



i A(<7_- — (T-| 



72 2e 



Therefore, the solution of (3.23) is 

n XiT ( C) = KJ!^±^± (l _ e -(— + )r) 



4 er_ — cr + 
i A((T_ — CT+) 



(i - e-nc - ' x ^--_ a+ \ i e-nc 



V2 2en V2 2e/Z 

Combining the solutions £(r) and f2A,r(C) with the expression for 7a,t(WXC)) 
pT9| we get 

1X,AW(0) = ^(e-^C)exp f-^j 2 g " +<7+ (l - e -<"— 



x exp 



72 2e/i 



\/2 2e/Z 
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To calculate £*(W(0) we use ( 3.6 1. Then for the first term one gets 
pS- 1 ( 7A , T (S^(())) 

= p e -,A( C +C)/V2 £ e ~Kb*-b)/e ( W (Qy«,'-b)/ t 



■ p exp 



ICI 2 ^- + ^+ 

4 er_ — cr_|_ 
1 ( A 



(1-e 



A 



., V/i (l-e-^X--(l-e-^)C) ) W(e~^0 



x exp 

If we put A = in this expression and substitute p for 1 — p, then we obtain 
the second term of (3.6). Together this yields the one-step evolution for the 
Weyl operator: 

\(\ 2 



Cl(W(0) = W(e-» T exp - 



jpcxp f 



4 a- - a. 



-(1-e 



(3.24) 



1 /A (i-^ r )C-ii-Oc))+i-p} • 



Note that operator £* is a convex combination (3.6) of quasi-free (completely 
positive) maps and hence, in general, it is not quasi-free itself. Evolution of 
the Weyl operator for the time t = tit (n-step evolution) follows from (3.24) 
by induction: 



(cn^(c)) = ^( e -^c)ex P 
1 ( x a 



4 a- - a. 



-(1-e 



— n((T_ — cr_[_ )r 1 



(3.25) 



x J] {pexp (~ Q(l - e -^) e -*^C - ~(1 - e-^e-^c)) + 1 - p}- 



>/2 \M M 
If <t_ > cr+j then Re(/^) < 0, which implies 

w* - lim W{e- nJiT C) = 1 



(3.26) 



see Appendix A. 3. To prove the convergence of the product (3.25) for the 
limit n — > oo, we denote 



^fc(C) :=p(exp 



1 (\ 



V2 Vm 



(1 - e-^ T )e- k ^C - »(1 - e-^ T )e-^ lT C ) - 1) 



A 



-flT\ -,-kfJ.T, 



(J- 



Since the product Y['kLo^~^^ lk (0) converges if and only if we establish conver- 
gence of the series J2kLa l^fc(OI' we nave to estimate the terms {|/ifc(C)l}fc>i- 
Note that 



-(a 



1/4 



)-r/2j e -k(a_-a + )r/2 



(3.27) 



Hence, (3.27) ensures the convergence of the series and the infinite product 
for <r_ — cr + > 0. 

Summarizing, we obtain that for the open cavity with parameters ct_ > 
<r_|_ > the characteristic functional (see Appendix A. 5) for the limiting state 
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w c,cr(') (1-50) exists and is given by 

ICI 2 



1 



4 er_ — (7+ 



n{ 

k=0 



/-lM^Qd-e^)e 



A 



CT^-tMT/ — — (1 



,-HT\ -kflT 



which is independent of the initial cavity state pc- 



(3.28) 
□ 



Remark 3.5. Notice that the evolution of the Weyl operator for the time 
t = tit can be written as a convex linear combination of quasi-free completely 
positive maps (see Appendix A. 4): 



(C) n (^(C)) = ^(e-"^C)exp 



4 cr_ 



+yj _ e -n(cr_-<T + )r 



n-1 

m=0 l<fci<...<fe m <n— 1 



N,, M>1 (Ml-e MT )(e 



— kifir 



-k m flT 



K-cc. 



Here c.c. stands for the complex conjugate of the first term. Then ( |3.28 ) and 
the last formula suggest that the limiting state ujc.a(-) is not quasi-free. 



The rest of the chapter is devoted to the proof of Theorem [L8j which is 
similar to the case — 0. The number of photons (1.55 1 for the time i = nr 
can be calculated using the adjoint operator 



N a (nr) = Tr c (b*b C n a { Pc )) = Tr c ((£* ) n {b*b) Pc ) , 



(3.29) 



where pc is the initial gauge-invariant state of the cavity. Note that by (3.6) 



and (3.7 1 we obtain for A = b*b 



C* a {b*b) = p £T 1 (e TX *.<' 5(6*6)) + (1 -p)e TL ^(b*b) 



(3.30) 



Lemma 3.6. The action of the adjoint operator C* on the photon number 
operator is explicitly given by 



Cl(b*b) 



(3.31) 



= e - {a —' J + )T b*b + p—e- (a — a +^(l - e» T )b* -p l -^e- {a — a+)T {l - e» T )b 



A 2 

p ri2 e 



>— <T +) T (l-e' iT )(l-e' iT ) + 



-(1-e- 



where p, = ie+ (er_ — o + )/2 



Proof. We start with the fi rst term in the right-hand side of (3.30). Since 
7a, t (^4) = e rL ">--°(A) (3.181, one can calculate j\ >T (S(b*b)) by taking into 
account (3.8). Then 



L* Xa ((b* - X/e)(b- A/e)) = iXb - i\b* - (cx_ - a + )(b* - X/e)(b - A/e) + cr+ 
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and 



Therefore, we obtain for mapping 7a,t(") the following system of differential 
equations: 

d^ T ((b* -X/e)(b-X/e) 
dr 

= -(<T- - (T + ) lx ,r((b* - X/e)(b - A/e)) + iA 7AiT (fc) - iA 7A , r (6*) + cr+ 

dj\,r{b) , . A((T_ - (7+) 

- -M7A,r(6) + - 



dr 
djx,r{b 
dr 



-fXj\, T (b* 



2e 

A(cr_ — er^ 



2e 



The solution of this system is 



7A,r(&*) = e^ r fo* + -(1 - + i *(1 - e"" T ) 

e 

3 -* T )-i^(l-e-* T ) , 
,«A 



7A,r(6) = e-^ T 6 + -(l 



(3.32) 
(3.33) 



7A,r((&* - A/e)(6 - A/e)) = e-^— a ^ T b*b + (— (1 - e^ 7 ") - -) e -^-- CT +H* 

/i < 

+ (- = (1 - e^ T ) - -) e -( CT -- ff +) T 6 + ^(1- e -(»— °4>) 



_ P -( ff --°'+)' r 



— 
cr_ — cr + 

A 2 (cr_ — fi+ ) sin er 



) + 



A 2 



5 -( cr -- CT +) T 



-(<r_-<r+)T/2 



(3.34) 



Making the shift transformation (2.1| of 7a,t((&* — ^/ e )(b— A/e)) and calcu- 
lating the second term in (3.30) by setting A = 0, we obtain (3.31 ). □ 



(3.321, (3.33) yield for k > 1 



If one plugs in (3.6) and in (3. 8), (3. 9) operators A — b* or = b, then 

(3.35) 



iX 



(C*) k (b*) = e- kp - T b* +p— (1 - e- kf * T ) 

A* 



(£* ) k {b) = e~ kfiT b - p— (1 - e- fe H 
A* 



(3.36) 



Since ££(•) | CT _ =(T+= o= £*(•), formulae (3.35), (3.36) coincide for er_ = cr + = 
with 12.191), (I2.20I) for the ideal cavity. 



Notice also that (|3_6| , (|3_7j and ([3^ , (^ imply 

=0. 



(3.37) 
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Proof, (of Theorem 1.8) To this end, we construct first a 4 x 4 matrix £* 
acting on the complex linear space spanned by the operator-valued vec- 
tors ( 6* 6,0,0 , ), (0, 6*,0, ), (0, 0, 6, 0), (0, 0, 0, 11)), according to the formulae 
|3.31 ), ( |3.35| , |3.36| and ( 3.37[ ), for k = 1. Then diagonalisation of £* al- 
lows to calculate powers (£* )™ an d to find explicit expressions for the n-step 
mapping 



(£*)"(6*6) = e- n{a -- a + ]T b*b + p—{e- nirT -- a ^ 7 

P 



>* (3.38) 



A* 
A 2 

H 2 



-(<7_-<T + )t/ 



e v— ^-(l- e ^)(l-e^) 



1 _ e -"(a'--o- + )r 
1 — e - ^- - ^) 1 " 



,2A 2 1 - e -n(<r--<r+)r 

p2 i_ e -("— (1 

9X 2 

~ (1 _ e -n(<r_-a + )r/2 



ImI 2 



cosner 



0"+ 



(1-e- 



Note that fl3.38[) reduces to ( |3.31| for n = 1. If r = 0, then (|3_38j) yields b*b for 
any n. By (3.38) we obtain for a gauge-invariant initial state the mean- value 

(3.39) 



of the photon number (1.55) in the open cavity at t = tit: 



n (h*h\\ - c ,-<°--°+)t 



N„(nr) :=Tr c (pcOCr(6*6)) = e 
+ p(l -P)ria(l - e" (CT - _ ct+)t/2 coser) 



iV<r(0) 
1 — e -( cr -- IT +)' r 



,2A 2 



P 



(1-e 



-n{(T - — a +)t / 2 



cosner) 



(1 



-n(cr_ — ct-)-)t 



Note that (3.39) coincides with (1.57) for t = nr. Finally, by virtue of (3.38) 
and (3.39) one gets for the u;*-limit (see Appendix A. 3) 



w* - lim (£*)"(6*6) 

n— >oo 

,2\ 2 1 - e -(°--- CT +W 2 coser , A 2 

= p(l — Pj-r—rK ; ? n 1- P(2p — l)r~P7 + 

y |/i| 2 l-e-( a -- a +> |m| 2 



(3.40) 



Therefore, lim^^ w c lr ff (6*6) = lim.^ Tr c (p c (£* ) n (6*6)) and ( |3.40 | yield 
the proof of (|1 .58b. □ 



Notice that in the limit of the ideal cavity: <j_ — » +0, ct_ > er + > 0, 
>+ o N a (t) — N(t), which is ( 1.28 ). For other limiting cases 



(3.39 1 gives lim CT __ 



see Remark 1 1.9 



Non-equilibrium states of a photon cavity 



25 



4. Energy and Entropy Variations 
4.1. Energy variation in the ideal cavity 



Since the time-dependent interaction in ( 1.3 ) is piecewise constant, the system 
is autonomous on each interval [(n— l)r, tit). Therefore, there is no variation 
of energy on this interval, but it may jump, when a new atom enters into 
the cavity. Note that although the total energy corresponding to the infinite 
system (1.5 1 is undefined its variation is well-defined [BJM1| -[BJM3 . 

The times when the n-th atom is actually traveling in the cavity are of 
the form t = n(t)r + with n(t) = n — 1 and v(t) € [0, t), see (1.15). 
To calculate the energy variation we compare two total energy expectations: 
for the moment t n — (n — l)r + v(t n ), when the n-th atom is present in the 
cavity, and for the moment £ n _i = (n — 2)r + v(t n -i), when the (n — l)-th 
atom was in the cavity. Then by ( |1.5[), ( 1.6 ) and by (1.16) we obtain for the 
total energy variation in the system (1.5) between two moments t n -\ and t n 
the following expression: 



4" (#(;„)) 



A£(i„,i„_i) .— lu s y±±\b n)) — uj s 
= Tr( e - iv ^ H "p s {{n-l)T)e- iv ^ H "H n ) 

= Tr(p s ((n-l)r)H n ) 

_ Ti(e- iTH "- 1 p s ((n-2)r)e iTH ^H n _ 1 ) 

= T^{T( n -i)T,o(PC ® PA)[H n - H n -x\) . 



(4.1) 



Here we used that the system ( 1.5 ), (1.6 ) is piecewise autonomous with H(t- 
0) = H(t), and the state ps{t) is u>*-time-continuous (Appendix A. 2). 



Recall that by duality (1.34) we have 



Tr (T(„_i) T)0 (pe ® PA)[H n - 
- Tr ( Pc ® PA T* {n _ l)Tfi {H n - H n ^)) 



Since (1.6) implies 



H n -H n _ 1 = A (b* + b) <8> (r?„ 



by ( |4Tj ), ( |42j ) and by [H v ,r)k] — we obtain 
&E{t n ,t n - X ) 

= Tr{ Pc ® p A T ( ;_ 1)T:0 (A (6* + 6) ® ® (ry„ - ?7n _ 1 )]} 
Lemma 4.1. For an?/ n > 1 one gefc: 

A 



(4.2) 



(4.3) 



(4.4) 



T* (b* ® 11) 



11 - 



-a 



fe=i 



T* r i0 (6 ® 11) = e - mTe b ® 11 - -(1 - e- ire ) ^ e -(n-fe)ir 611 , 



(4.5) 



(4.6) 



fe=i 
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Proof. Let us define 



B* k (T) :=T* k (b*®t) 



(b* ® l)e 



-itH, 



, k> 1 



(4.7) 



Then by (1.61 and ( 1.32[ ), (1.331 the operator (4.7) is solution of equation 

d s B* k (s)=i[H k ,B* k (s)]=ieBl(s)+\l®ri k , B* k (0) = 6* <8> 11 , 
which has the following explicit form: 



B* k (T) = e tTe (b* (8)1)- -{l-e lTe )l® Vk . 
e 



Similarly one obtains 



Note that by iteration of (4.8) for k — 1, 2 one gets 

T 2 * T)0 (&* ® 1) = e 2lTe (fo* <g> 1) - -(1 - e lTe )e lTt t® m 



(4.8) 
(4.9) 
(4.10) 



A 



(l-e lT£ )10 77 2 



Proceeding with iteration of (4.10) we obtain (4.5). Using (4.9) one proves in 
a similar way (4.6). □ 

Recall that in the present paper we suppose that atomic beam is homo- 
geneous \\.2\ , i.e. p = Tr{pc®PA(^Vn)} is independent of n probability that 
the n-th atom is in the excited state, and that the atomic beam is Bernoulli 
(Remark [TTT]) : 

Ti A {p A (r] ni T] n2 )} = S nun2 p+ (1 - S nun2 )p 2 . (4.11) 



Let the initial cavity state pc be gauge-invariant. Then by Lemma |4.1| and 
( 4.11[ ) one obtains: 

Tr{ Pc ® PA T ( * n _ 1)T)0 (A (b* +b)(g> l)[l<g> {r, n - Vn ^)}} 



(4.12) 



-(1 - e ire ) e {n - k - 1)irer ftAPAVk(Vn ~ Vn-l) 



(1 - e 



fe=i 

n-l 



, — (n— k— l)ire 



2A 2 

p(l — p) (1 — cos re) 

£ 



Hence, formulae (4.4), (4.12) prove for the energy variation on the interval 
(t n _ i,t n ) the following statement. 



Theorem 4.2. For the case of the ideal cavity the total energy variation (4-1) 
between two moments t n —\ and t n , where n > 1 is 



2A 2 

A£(t n ,t n -i) = p(l -p) (1 - cos re) , 

e 



(4.13) 
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i.e. for the variation between to = v{t$) and tn > to we obtain: 



2A 2 



A£{t n ,t ) =J2AS{t k ,t k - 1 ) = (n - 1) p{l - p) (1- cos re) . (4.14) 



fc=i 



Remark 4.3. The total energy variation (4.1 ), when the n-th atom is traveling 



through the cavity between the moments t' = (n — l)r and t" = nr — 0, can 
be written as 

A£(t", if) = u^(H n ) - 4 n " 1)T (^n) ■ (4-15) 



Here again we used that for t € [(n— l)r, nr) the Hamiltonian ( 1.5 ) is piece- 
wise constant and that it has the form (1.6), as well as that the state w^(-) is 
ui*-continuous in time (Appendix A. 2). Since the system ( 1.5 1 is autonomous 
on the interval [(n — l)r, nr), one obtains A£(nr — 0, (n — l)r) = 0. Then 



(4.15) implies that on this interval (in contrast to (4.4)) the variation of the 



interaction-energy completely compensates the energy variation due to the 
photon number pumping: 



LUI T (\ (b* +b)® 7] n ) - 4" 1K ( A ( b * + b )® Vn-l) 

= -[u)f{eb*b® 1) -Lu ( s~ 1)T {tb*b® 1)] . 



(4.16) 



Note that similar to (4.12 ) one can check this identity explicitly using Lemma 



4.1 and (4.11) applied to the left-hand side of (4.161 



4.2. Energy variation in the open cavity 

Although for the open cavity the time-dependent generator (1.37) is still 
piecewise constant ( |1.38 1, the cavity energy is continuously varying between 
the moments {t = kr}k>o (when the interaction may to jump ( 1.5 1) because 
of the leaking/injection of photons. 

Therefore, as above we first concentrate on the elementary variation of 
the total energy, when the n-th atom is traveling through the cavity between 
the moments t! = (n — 1)t and t" = nr — 0: 

AS a (t",t') := u^ a {H n ) - 4"; 1)r ( J ff„) . (4.17) 
Here again we used two facts: (1) for t E [(n — 1)t, nr) the Hamiltonian 



(1.5) of the form (1.6) is piecewise constant; (2) the state 



(1.42) is 



time-continuous (Appendix A. 2). 

By virtue of ( |1.42[ ) and of ( 1.43 1 for operator A = H n ( 1.6 1, we see that 
the problem ( |4.17[ ) reduces to calculation of the following expectations: 

e 4 T a(&*& ® A) and A w|y (&* + b) ® r, k ) , k, s > 1 . (4.18) 



(3.391 



The first expectation in (4.18) is known due to (1.46) and Theorem 1.8 



.At 



c u, 5iff (&*6® 1) = e N a (kr ) . (4.19) 

To calculate the second expectation in ( |4.18 1 we use (1.42) for the operator 
A = ((b* + b) <E> 1)(]1 £g) rjk) and the representation (1.43) for initial gauge- 
invariant state pc and for homogeneous atoms state pa- 
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Lemma 4.4. Let ct_ > <r + > 0. Then for the mappings {{T^ r Q )*} n >o, see 



1.43), one obtains: 



(T^ r0 )*(6® 1) =e- n » T b® 1- -(l-e-" r )V e- {n - k) ^t® m , (4.20) 



k = l 



and {T° )*(b* ® 1) = ((T- )*(6® !))• 



Proof. We prove this lemma by induction. Suppose that formula (4.20 ) is true 
for (T^ t0 )*. Then we show that it is also valid for (Tf +1 v T )*. By virtue of 

(1.431 and by (3. 4), (3. 11 1 for the action of operator e CT >™, one gets 
(T n CT +i r(&®fl) = (T„T(e^.» +I (6®ll)) 



(^)*(^ 1 (e7^" +1 (5 n+1 (6®l)))) 

e 



(^)*(^7+i(e TZ; - +1 ((6 - ~) ® r?„+i + 6 ® (/ - fj»+i)) )) , 



where we used (2.2) in the last line. By (3.131,(3.14) and (3. 18), (3. 33) com- 
bined with the shift S~l t (2.8), we obtain 



S-^(e r£ -»^((b- -) ® Ifc+l)) = %((7A,r(6) - ") ® 

= (e^ T 6--(l-e^ T ))®?7n+i , 
A 4 

^(e^+^fe® (7 - ifen))) = S-^(7o,r(6) ® (I- Ib+l)) 

= e-^6® (J-tJn+l) ■ 
Consequently, 

= (T„ ff )*((e-^6 - -(1 - e^)) ® 7? n+1 + e~» T b ® (1 - ?7 n+1 )) 
= (T„ CT )* « 1 1 - -(1 - e-^ T )]l ® ry n+1 ) 

n+l 

= e - {n+1) ^ T b ® 1 - — (1 - e^ T ) V e -(n-*+i)/^j ^ %j 



fc=l 



which proves the lemma. 



□ 



Recall that p = it + (<t_ — o + )/2. Hence, in the limit cr_ — > +0 one 
recovers from this Lemma formulae (4.5) and (4.6) for the ideal cavity. 

Corollary 4.5. Let initial cavity state pc be gauge-invariant state for homo- 
geneous state pj^ of the atomic beam. Then with help of (4.11) and (4.20) 
one obtains the interaction energy expectations (4.18) corresponding to the 
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difference (4.17): 



\u,^ a ((b*+b)®r ln ) = 
2A 2 e 



(4.21) 



p(l - p)(l - e-^-^+^Z 2 COSer) + p 2 (l - e -n(<r—<r+)r/2 CQf , ^ 



IMI 

A 2 (<t--<t + ) 
H 2 



p(l - p)e- (CT -- <T + )T/2 siner + p 2 e -^--<r + )r/2 ginner 



A 4"; 1)T ((6* + 6) ® 7?„) = -y^- P 2 (l - e -(«-i)("— C0S („ _ i) er ) 

I Ml 



A 2 (CT- - CT+) 



iMl 



p2 e -(n-l)(<r_-a + )r/2 sin ( n _ ^ 



(4.22) 



Corollary 4.6. Taking into account Theorem 1.8 and (4.21 ), (4.22 ) we get for 
the elementary variation of the total energy (4.17) 

A£ a (nr - 0, (n - l)r) = e (iV ff (nr) - N a ((n - l)r)) (4.23) 

+ a k; ct ((6* + 6) ® - 4"; 1)r ((&* + 6) ® ^)) 



- AU0))(1 - e -(»--<^+)T) e -(n-l)(«r— a+)7 



2A 2 e 



>_-<t + )t/2 



coser)(l — e 



-(n-l)(<T_-(T + )r 



+ p(l-p) 

\2 



A 2 ((7_ — (7 



+Z e -( CT _- CT+ )r/2 sineT 



P 



2 A 2 (cr_ -<7+) 



Ml 



e -„( CT _-<r + )r/2 sinn£T _ e -(n-l)(«r_-«r + )r/2 



sin(n — l)er 



Note that in the limit of the ideal cavity: c+ — >• and cr_ — > 0, one finds 
for total energy variation (4.23): Af CT (nr — 0, (n— l)r) = 0, which corresponds 
to the autonomous case, see Remark |4.3| Whereas for <r_ > er+ > the 



external pumping due to a + /(er_ — cr + ) > is in competition with the energy 
leaking, see the second term in the right-hand side of (4.23). The limit of the 
energy increment when n — > oo is 



lim A£(j (tit — 0, (n- l)r) 



(4.24) 



A 



2 r 



p(l-p)~ -2e(l -e"<>— ^^coser) + (tr_ - cr+Je-^--^/ 2 si 



sin er 



IMI 



Remark 4.7. To consider the impact when the n-th atom enters the cavity we 
study the total energy variation on the extended interval {(n— l)r— 0, nr- 0). 
Then 

A£ ff (nr - 0, (n - l)r - 0) = (^""(fl,) - 4"; 1)T (ff„)) (4.25) 

+(4:; 1)T (^)-4:; 1)T_0 (^-i)) ; 
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where the second difference A£ a ((n — l)r, (n — 1)t — 0) 



(n-lV 



(Hn) 



(n-l)r-O 

O ,(T 



(i/„_i) corresponds to the energy variation (jump), when the n-th 



atom enters the cavity and the (n — l)-th atom leaves it. 

To calculate A£ a ((n — l)r, (n— 1)t — 0) note that by the time continuity 
of the state 

A£ a ((n - l)r, (n - 1)t - 0) - Tr(p 5 ((n - l)r)tf n ) - Tr(p 5 ((n - l)r)# n _i) 
= Tr (e^-- 1 ... e^Goc ® ^))(i/„ - 

= Tr (l£_ 1)Ti0 (pc ® AU)(#« - #»-l)) 

= Tr ( Pc ® ^(T^_ 1)Ti0 )*(A(&* + 6) ® - ?7«-i)) 

= Tr ( Pc ® ^(r ( CT n _ 1)T!0 )*(A(6* + 6) ® (% - ?7n-i)) 
= Tr{ Pc ® P^(T ( CT „_ 1)T:0 )*(A(6* + 6) ® 1)[1® (r/„ - 



where T£L nT0 = e rL — ... e Ti " 1 is defined by (1.40). 

If the initial cavity state is gauge-invariant, then (4.201 yields for the 
energy jump at the moment t = (n — 1)t: 



A£ CT ((n-l)r, (n-l)r-Q) = 

Tr{p c ® Px(r ( CT „_ 1)Ti0 )*(A(6* + 6) ® l)[t ® (??„ - ??n-i)]} 



71-1 



(1 - e"^) ^ e- ( "-' £ - 1) ^Tiu(^%( ?7 „ - Tj n — i ) ) 
fc=i 
n-l 

(1 - e"^) ^ e-^-'-^^^Cp^^CTM - Vn-i))- 



(4.26) 



fc=i 



Taking into account the Bernoulli property (4.11) we obtain from (4.26) 

A£ a ((n - 1)t, (n - l)r - 0) = — (1 - e"^)p(l - p) - 5(1 - e"^)p(l -p) 
2A 2 e, 



: p(l - p)tV(! - e-( CT — CT +) T / 2 coser) 
A 2 (ct_ - a 



p(l -p)- 



±Z e -(«r— a +)r/2 sineT _ 



IMI 



(4.27) 



Notice again that for <r_ 



-0 one obtains from (4.27) the one-step 



energy variation for the ideal cavity (4.12). 

Summarising (4.23) and (4.27), we obtain the energy increment (4.25) 
which is due to impact of the open cavity effects (4.23) and to the atomic 
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beam pumping (4.27): 
A£ a (riT-0, (n-l)r-O) 



(4.28) 



cr_ — a j 



7V CT (0))(1 - e -(»— o-+)r) e -(n-l)(«r_-o-+)T 
(<x_-<x + )t/2 coser ^-(„-l)( (7 _- CT+ ) T 



, A 2 (g_ -<7 4 

IH 2 



e -„(a_-<r + )r/2 sinn£T _ e -(„-l)(«r_-«r + )r/2 



sin(n — l)er 



Theorem 4.8. -By virtue of {4-28) the total energy variation between initial 
state at the moment to := — 0, when the cavity is empty, and the moment 
t n := nr — 0, just before the n-th atom is ready to leave the cavity, is 



A£ a (t n ,t ) = J2 Af CT (fcr - 0, (k - 1)t - 0) 



(4.29) 



fc=i 



= e(^ T± 7V CT (0))(l-e-" (,T -~ CT + )T ) 

<7_ — (T_|_ 

+ p (i _ p )±li(i _ e -(*— <* +W 2 cos e r) " 



IMP 



1 — e -( cr -- CT +)' r 



_ffere -/V^O) = lo^ 1 a (b*b ® 1) is t/ie initial number of photons in the cavity. 



Remark 4.9. Note that the total energy variation (4.29) is due to evolution 
of the photon number in the open cavity (3.39) and the variation of the 
interaction energy (4.21), that give 

A£ a (t n ,t ) = euj^ a {b*b ® 1) + A w£,((&* + b) ® 77„) (4.30) 

- ec4 )O .(&*&® 1) . 

For <r_ — cr + > it is uniformly bounded from above 

a+ 2A 2 e p(l - p) 



A£ a (t n ,t ) < e 



cr_ — cr + |/i| 2 1 — e 
2 A 2 ((7_ -CT+) 



-(<7_- CT+ )T/2 



(4.31) 



I 12 

H 2 



The lower bound of (4.29) is also evident. It strongly depends on the initial 
condition ^(0) and can be negative. 

The long-time asymptotic of ( |4.29 1, or (4.30), is 



A£ a := lim A£ a (t n ,t ) = e( 

n— » — 

+p(i-p) 



<7_ — CT-| 

2A 2 e 1 - e -( <7 -- <7 +) T / 2 coser 



- JV ff (0)) (4.32) 



1/4 



1 — e -( CT --°~+) T 
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From (4.32 1 one gets that in the open cavity with cr_ — cr + > the asymptotic 



of the total-energy variation is bounded from above and from below 
(J+ „ 2A 2 e p(l -p) 



Af ff < e ( 
A£ ff > e 



<7_ — (T_| 
' <7_ — (7 4 



-iV CT (0)) 



| M | 2 l 
2A 2 e 



e -(cr_-<T + )r/2 
P(l -p) 



|/X| 2 1 + e -(o---«r+)r/2 



For the short-time regime nr <C 1 one gets for (4.29) 



A£ CT (i n ,i ) = nT€a + - nr(a- - cr + )iV (T (0) 



(4.33) 
(4.34) 

(4.35) 



+ nrp(l — p)-: — ry(l — coser) 



(T_ — (Tj 



+ nrp' 



W 

A 2 (fi_ -g + )e 



+ G((nr) 2 ) 



1 — e -( CT -- cr +) 

2\ 



i.e. a linear asymptotic behaviour. 

Another asymptotics one finds for the small difference between leaking 
and external pumping: er_ — <r + — > 0. Then (4.29) yields linear behaviour 

(4.36) 



2A 2 



A£ CT (t„, to) = nTe cr_|_ + np(l — p) (1 — cos er) , 

which is a growing of the total energy due to the both external and atomic 
beam pumping. Note that in this limit the energy variation A£ a (t n ,to) is not 
bounded from above (4.31), (4.33). 



This case coincides with result for the ideal cavity (Theorem 4.2) when 
the rate of environmental pumping a + = 0. 

4.3. Entropy production in the ideal cavity. 

One of the central quantities to study in the non-equilibrium statistical me- 
chanics is the entropy production (or the entropy production rate). We refer 
to the series of papers |BJMl] - |BJM3j by Bruneau, Joye, and Merkli, for a 
detailed discussion of this quantity in the context of open quantum systems 
with repeated interactions, and we shall adopt definitions of these authors. 

Let p and po be two normal states on the algebra 21(H). We define the 
relative entropy Ent(p|p ) of the state p with respect to a "reference" state 
Po by 

Ent(p|p ) :=Tr(plnp-plnp ) >0. (4.37) 

The non- negativity follows from the Jensen inequality: Tr(p In B) < In Tr(pB) , 
applied to observable B := po/p. 

Here we calculate the entropy production: 



AS(t) := Ent(p s (t)|^ e/ ) - Ent(p s (0)\p r s ef ) 



(4.38) 



for the ideal cavity: er_ = cr+ = 0, with dynamics ps '■ t i-> ps(t) (1-16), 
and for a reference state p r ^ <X> pa- To make a contact with thermodynam- 
ics, we suppose that all atoms of the beam are in the Gibbs state with the 
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temperature 1//3. Formally this can be written as 



-PH An 



see (1.5). Since 



P r s f = p r c S ® PA = {Pc f ®m® PA) , 



and since for the unitary dynamics (1.9) 



Tr { Ps (t)\n p s (*)} = Tr{p s (0)lnp s (0)}, 



the relative entropy (4.38) is 

AS(t):=Tr e {[p^-pP]lnp r e ef } 



(4.39) 

(4.40) 
(4.41) 
(4.42) 



J2^{[Ps(0) - p s (n T )](t® H Ak )} , 



fc=i 



for t = tit + v 1 see (IL15J). Here is defined by \l.22) . 

Remark 4.10. For any Hamiltonian H n that acts non-trivially on He <8> H An 
and for any Hamiltonian H Ak acting on H Ak one gets [H n , H Ak ] = if n ^ k. 
Note that in our model ( 1.21 1 implies [H n> H Ak ] = for any n,k. 



Therefore, by virtue of (1.22) and (1.23) we obtain 



Tr{e- lTH "...e- lTHk + 1 p s (kT)e lTHk + 1 



Tr{p s (nr) (l®H Ak )} 

-irH n ^ — irHk 

= Tr{p 5 ((fcr) (l®H Ak )} ■ 
Then by the same arguments one gets also that 
Tr{ Ps (Q) (l®H Ak )} 



(4.43) 



JtH„ 



(4.44) 



Tr{e- tT±lk - 1 ...e- lTMl p s (0)e 
Tr{p s ((k-l)r) (l®H Ak )} 



irHk-i 



(*®H Ak )} 



In the case of our model (see Remark 4.10) H Ak = e lTHk (H Ak )e %rHh . 
Therefore the last formula gets the form 

Tr{ P5 (0) (l®H Ak )} = Tv{p s (kr) {l®H Ak ). (4.45) 

Equations (4.43) and ( 4.45[ ) shows that the second term in the entropy pro- 
duction (4.42) vanishes. 

If we suppose that the reference state is the Gibbs state (1.29) at tem- 
perature l/flc , then by (1.26) one gets 

AS(t)=Trc{[Pc° ) -Pc n) }tePc f } (4-46) 
= Tr c {[p( 0) - p W](-(3 e eb*b)} 
= Pce(N(t)-N(0)) , 



where N(t) is the mean photon number defined by (1.26). 
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Let us define by A£ c (t) = e(N(t) — N(0)) the energy variation of the 
cavity due to the photon number evolution. Then (4.461 expresses the 2nd 
Law of Thermodynamics 

AS(t) =/3 c A£ c {t) , (4.47) 



for the pumping by atomic beam. Note that relation (4.471 does not depend 
on the initial cavity state p£ 



(0) 



Remark 4.11. In general, when [H n , H An ] ^ 0, the combination of (4.42 1 with 
(4.431 and (4.441 yield for the entropy production at the moment t = nr-\-v 
the expression: 

AS(t): = Tt c {\pf -/4 n) ]m/£ e/ } 

n 
k=l 



(4.48) 



The last term in (4.481 can be rewritten into the standard form [BJMlj - 
[BJM3], if one uses the identities: 

Tr{ Ps (kT)(t(g>H Ak )} = Tr{e TLk ...e TL i( Pc ® PA ) (1® H Ak )} 



fc-1 

Tr{e-' T ^ (e TLk -\ . . e Til [p c ®(g) p n \) ® p k e lrH " (1 ® H Ak )}[l ® (g) p n 

n—1 m>k 
JrH k /-n u . \-irH k 



and 



Tr{p s ((fc-l)r) (l®^J} = Tr{4 fe X) ® pfc-EU J . 
For i = tit + this gives the formula for the entropy production in the ideal 
cavity: 

AS(t) = Tr c {[p( 0) -p^]lnp^} (4.49) 

n 

+ /3 ^Tr Wc0 ^ fe {(p^- 1) ®p,)[e i ^(l®^Je- i ^-]l®^ fc ]} . 
fc=i 

One expects that for an open cavity the gain and loss of photons would 
result in a corresponding additional flux of entropy due to the quantum 
Markov evolution. It is not yet completely clear (see e.g. jFaRe] ) how to de- 
fine this entropy production correctly and whether contributions of these two 
processes are independent. Therefore, we leave the analysis of the entropy 
flux for the open cavity to be considered in the future. 



5. Concluding Remarks 

From Theorem |1.2| we learn that in the ideal cavity the photon number in- 
creases linearly in time up to bounded oscillations, (see Figure[TJ green curve) 
The rate of the growth is non-zero for any value of the probability of excited 
atoms in the beam except p = and p = 1. For large time it is independent 
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Figure 1 . In blue, we have plotted the photon number in an 
open cavity with parameters cr_ = 0.003, a + = 0, e = 0.5, 

1 A 2 

r = 0.5, p = |, ]T^T2 = 1 and vanishing initial number of 
photons N a (0) = 0. The green curve is the mean- value of 
photons in the ideal cavity for the same parameters and ini- 
tial condition N(0) = 0, except er_ = 0. The expressions for 



these quantities are given by (1.281 for the ideal cavity and 



by (3.39) for the open cavity. The red line is the asymptotic 



value (1.58) of the mean photon number in the open cavity. 



of the initial state pc- The rate of the linear growth of the mean- value of 
photons N(t) with respect to the time t = nr can be seen from (1.28) since 



N(t) N(0) , . 2A 2 , . p 2 2A 2 , , , 

— ^ = — ^ + p i -p) — =■ 1 -coser + L — —(l-cosner) . 5.1 
tit nr re z nr e l 

Hence, for JV(0) = one gets linear growth modulo bounded oscillations: 



Nit) , . 2A 2 
_Li=p 1-p _(1 



p 2 2\ 2 

coser)H ^-(1 — cosner) . (5.2) 



This linear growth of the photon number is observed in experiments with one- 
atom masers for the high-quality resonators (nearly ideal cavities) |M WMj . 

As can be clearly seen in Figure [T] (blue curve) , the mean photon num- 
ber N a (t) in the open cavity ( 3.39[ ) is also initially increasing linearly up to 
bounded oscillations, but then it stabilizes for large times. Similarly to (5.1) 



the rate of the growth in the open cavity with respect to the time t = nr can 
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be deduced from 

N^nr) =c - n(a _- tr+)T N tr (0) 

TIT TIT 

o\2 i _ n(o-_— cr+)"7 

+ P(1 -P)r7a(l -e"( CT — CT +W 2 cos «-) - 



(5.3) 



|^| 2V 'nr(l -e-C "--^) 1 ") 

_ e -n(— + )r/2 CQsneT) + °+ L_g 

nr cr_ — cr_|_ nr 



which for <r+ = and then cr_ = gives ( |5.1[ ). 

Let iV^O) = and <r+ = 0. Since /i = (cr_ — o + )/2 + ie, the short- 



time behaviour of (5.3 1 for n(er_ — cr+)r < 1 is linear (modulo bounded 
oscillations) : 

NJnr) ,„ ,2A 2 . p 2 2A 2 

— ^ i =p(l - p ) — (1 - coser) + cosner)+ (5.4) 

TIT Te z TIT e z 

+ <D((n(a_ - <j + )t) 2 ) , 



and asymptotically it is close to (5.2 1. This is clearly visible in the Figure [I] 
for n < 100. 

There are several generalisations of the beam-cavity problem considered 
in this paper that could be handled by suitable modifications of our methods. 
First, there is the detuned case, when the distance between atoms is greater 
than the length of the cavity, d > I. This is a situation where there is still at 
most one atom interacting with the cavity at any given time, but there are 
time intervals without an atom present. The modifications needed to analyze 
this situation are straightforward. A more interesting generalization would 
be to consider random interatomic distances d > I, when again still there is 
at most one atom in the cavity, but they arrive randomly. 

Due to properties of the apparatus that produces the atom beam, one 
may expect short range correlations in the chain of atoms. Such a correlated 
atomic beam can be described by a classical Markov chain or, more generally, 
by a so-called Finitely Correlated State [FN Wj . Calculating the asymptotic 
behavior of the cavity in this situation will be a bit more complicated but 
should still be doable. 

We would like also to mention the following two other problems. The 
first is to consider a cavity interaction with 'soft atoms'. One of the possible 
interactions is of the form b* (8) a~ + b <8) a + , i.e., interaction of the Jaynes- 



Cummings type. Here in contrast to (1.3) the atomic operators a := (a x ± 



ia v )/2 are off-diagonal, constructed with the Pauli matrices o~ x and a v . For 
p < 1/2 and for an ideal cavity this problem has been studied in |BPi| . It 
was shown that the Jaynes-Cummings interaction allows the cavity state to 
converge to some explicit thermal state. So, the pumping is saturated and 
the photon number expectation is well-defined and bounded. 

The second problem is to calculate the entropy production for the 
Kossakowski-Lindblad open cavity and to determine its relation to the en- 
ergy flux in the frame work of our model. For a discussion of the entropy 
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production problem for the quantum Markov semigroups see, e.g., the recent 
paper [FaRe] . 
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Appendix: Open One-Mode Photon Cavity 

Here we collect some remarks and recall certain statements concerning the 
quantum theory of open systems, see e.g. [Davll lAlFal IAJPII] . To this end 
we consider example of the open one-mode photon cavity C. This soluble 
example is useful for our discussion of the leaky cavity pumped by atomic 
beam starting on Section [3] 

A.l Markovian Master Equation. We treat the cavity: He = eb*b, interact- 
ing with external reservoir 1Z in the framework of the Markovian approach 
[AJPIIl IAJPIII] . Then 1Z is a source of a leaking, which decreases the cav- 
ity energy with the rate <7_, or/and of a pumping with the rate a+. The 
corresponding Markovian master equation for evolution of the normal cavity 
states with trace-class density matrices pcif) € £i(%c) is extension of 

the Hamiltonian dynamics by Kossakowski-Lindblad damping (leaking) and 
pumping terms [ATI lAlFaj : 

j t Pc(t) = L c ,*{pc(t)) := -i[H c ,Pc(t)] (5.5) 

+ ±a-([bp c (t),b*} + [b,p c (t)b*}) + ±a + ([b*p c (t),b] + [b*, Pc (t)b]) ■ 

Here pc(t) € dom(L^ CT ) and parameters er_,er + > 0. 

Note that evolution for p c := p c (0) G t\{Hc)- 

4,„ : PC pc(t) , (5.6) 

is the case of trace-norm continuous semigroup: C f c a :— e tLc - a on the Banach 
space £i(%c), with unbounded generator Lc,cn see e.g. |Za| . Ch.2.4. 



It is known that this case of quantum Markovian dynamics (5.6) needs 
a special care, see e.g. |Davl| . [ST1 IChFa] . but for the open one-mode photon 
cavity C all necessary properties can be checked explicitly. 

Denote by Lc ;(T the non-Hamiltonian part of unbounded generator ( 5.5 ). 
Then for any p e <£i(Hc) one has: 

LcAp) = E °«{V a pV* - \(V:V aP + P V:V a )} , V X = b , V t = b* . (5.7) 

By virtue of the trace cyclicity this canonical form of the Kossakowski- 
Lindblad generator: 

L C , CT := -i[H c ,-]+Lc,a , (5-8) 



ensures the trace-preserving property of dynamics (5.5) 

d 



d Tr c p c (t) = 0. (5.9) 



Note that by virtue of (5.5) the Markovian dynamics (5.6) is also unity- 
preserving: C} c a (T) = 1 for t > 0. 

To check another important property: the (complete) positivity of the 
trace-norm continuous semigroup {C c CT }t>o on the space £i(Hc), let us 
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present its generator (5. 



as L, 



:= — r, where 



$0):= ^V aP V* , ct q >0 , (5.10) 

a=i,t 

r(p):=Vp + py* with *:=*ffc + i a <*V a V* . (5.11) 

a=4->t 

First we reduce our analysis to positivity and we postpone the question con- 
cerning complete positivity to the end of this section and to the Heisenberg 
picture of quantum dynamics (5.6 1, see subsection A. 2. 



To see that dynamical semigroup {e**} t >o with generator (5.10) enjoy 
the property of the positivity, notice that by (5.10) one gets for the trace- 

" K /<:<: 



continuous maps p p®(t) := 
d 



dt 



^ CT a V a p<s> (t) V* 



(5.12) 



Let p e dom($) C <£\{H C ) and p > 0. Then ( |5.10| ) implies that > 

and that equation ( 5.12[ ) is positivity-preserving. This yields positivity of the 
solution p<s>(t) for t > 0, if /?$(i = 0) = p. 

To see that semigroup {e~ tr }i>o is also a family of positive mapping 
on Ci('Hc) note that (5.11) yields 

d 



dt 



Pr (t) = -T(p r (t)) = -(*pr(t)+p r (*)**) 



d_ 
dt 



e- t9 pe- t9 



(5.13) 



Then the mapping e tr : p *-> e **pe ** is po sitive. For p £ Ci('Hc) we 
denote by pr(t) 
for p > 0. 



e (p) the solution of (5.13). This operator is positive 



e t *(e^ r 



By virtue of (5.12) and (5.13) the composition of two maps: F(t) : p H> 
(p)), is a positive trace-norm continuous mapping. Note that the 
powers of F(t) are also positive maps. Then it is also true for dynamical 
semigroup with generator (5.8), since by the Trotter product formula one 
gets 



4,a = e 



the 



\\ — lim 



e t*/n e -*r/j 



(5.14) 



in the trace-norm topology |Zaj . This remark ensures, in particular, the £i- 
continuity of the limit (5.14). Note also that by (5.5) and (5.6) the limit 
(5.14) is unity-preserving (or Markov) semigroup. 

Recall that essential in the concept of open systems is a coupling of 
some small sub-system with "environment", which is a certain large (even 
infinite) system. The mathematical description of this concept involves a 
tensor product of corresponding configurations spaces and states or algebras 
of observables. Then the positivity-preserving evolution due to the master 
equation for a coupled system, must be robust for forming the tensor products 
[Davl]. Since tensor product of two positive maps might fail to be positive 
( [AlFaj . Ch.8.3), the evolution (5.6) has to verify a stronger condition than 
positivity established in (5.14). 
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We recall now the complete positivity property and constrains that 



it implies on dynamical semigroup (5.6| and on its generator [Davli lAlFal 
IA.IIMH. 

Let T : — > 21 1 - 2 ' be a positive linear map between two C*-algebras, 
i.e. T(A) > 0, where A G 2t (1) and A > 0. For k = 1, 2 and for n G N, let 

M n {% {k) ) ~ 2l (fe) ® M(C n ) , (5.15) 

be algebra ofuxn matrices with entries in 2l( fc ). Each of A4 n (f2L^) is also a 
C*-algebra. If we denote by Id„ the identity matrix from A4(C n ), then 

T n := T ® Id„ : X„(2l (1) ) -> M„(2l (2) ) , (5.16) 

defines a linear map by acting with T on each of the matrix element of 
the operator-valued matrix A n G ■M„(2l < - 1 - ) ). The positive map T is called 
n-positive (respectively completely positive), if operator T n is positive (re- 
spectively (5.16) are positive for all n > 1). For n = 1 it obviously reduces to 
the positive map. 

Example 5.1. A simple example shows that this property is quite non-trivial. 
(For more of them we refer to |AlFa] .) Let 2l^ =1 ' 2) = M(C 2 ) be C*-algebra 
of square complex matrices. Then the map of matrices A^(C 2 ) to the adjoint, 
T a dj ■ A — > A* is obviously positive. Denote by {Eij e A / J(C 2 )}i J= i i 2 the 
set of matrices with 1 in the ij-th entry and zeros elsewhere, i.e. Id 2 = 



algebra (15.15) 



=1,2 



To verify whether the map T a dj is 2-positive we consider the 



(5.17) 

E* and 



M 2 (M(C 2 )) 

and the element E 
E 2 = 2E, it is positive E > 0. On the other hand by definition (5.16) one 
gets 



M(C 2 )(g>M{C 2 



X(C 4 ) , 
C 4 ). Since E 



Tadj, 2 '■= T a dj ® W2 : E 



Tadj (En) 
Tadj (E21 ) 



T a dj (El2) 
T a dj(E22) 



(5-18) 

The matrix T a dj,2(E) in (5.18) is not positive, since its spectrum contains 
(—1). Therefore, the map T aa j is not 2-positive either. 

One of the important corollary imposed by demand that quantum Mar- 
kovian dynamics on the Banach space £i("H): 



P^ P {t) , p(0)= P e£i(U) , 



(5.19) 



must be completely positive, is certain restrictions on the form of the genera- 
tor L, cf. (5.6) and (5.8). For the case of bounded generator, when semigroup 



(5.19) is continuous in the operator-norm topology of mappings on £i(H), 
this is just a celebrated Kossakowski-Lindblad result saying that the most 
general form of L is 

L := -i [H, •] + L , (5.20) 
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where non-Hamiltonian part can be presented as 
1 



Hp) 



2^ 



{[V aP ,V*] + [V a ,pV*]) 



(5.21) 



-(V*V aP + P V*V a )} 



see, e.g., [ AlFaj , Ch.8. Note that t he choice of bounded operators H = H* 
and {V Q } a in representation (5.21) is not unique. 

The co ntact b etwee n (5.20), (5.21) and representations (|5.7[ >,( 5.8 ) as 
well as with (5.10), (5.11) related to semigroups {e**} t > , {e~ 4i }t>a follows 



through verbatim. Then taking into account that the Stinespring theorem 
([AlFa], Ch.8) implies the complete positivity of these two semigroups and 
using the Trotter product formula one can check the complete positivity of 
( 5T9 I- 

The present results for unbounded H and/or L are more poor. Certain 
classes of strongly continuous on contraction semigroups (5.19) have 

been studied in [Davlj and [EvLeli IEvLe2] under condition that operator 
J2 a V<y.Va is generator of a strongly continuous contraction semigroup on H. 
For further developments see, e.g., [Dav2 [ ISil IHoll IChFaj . 

The Kossakowski-Lindblad representation (5.20), (5.21) explains our 
choice of the right-hand side in the Markovian Master Equation (5.5), but 
appeals for a concrete verification of the complete positivity of quantum 
Markovian dynamics (5.6 ). Similar to other known cases of the Kossakowski- 
Lindblad type generators |FrVe| , this property of L\> a follows directly from 
explicit calculations. 

A. 2 Dual Dynamical Map (Heisenberg Picture). The equivalent (and often 
more convenient) is the abstract version of the reduced Markovian dynam- 
ics A.l on the algebra of observables 2l(7ic), i- e - the quantum dynamical 
semigroup ££.* CT := CT )* in the dual Heisenberg picture | AlFa| lA"JPl] . 

Recall that in a general setting the C*-dynamical system is a pair (21, r') 
with two properties. First, 21 is a unital C*-algebra, i.e. 1 € 21. Second, r* is a 
strongly continuous one-parameter *-automorphism of 21, i.e. for any A E 21 
the map: 1 1-> r*(A), is continuous in the norm topology of 21. 

Since quantum states belong to trace-class <Zi(H), which is not a uni- 
tal C*-algebra, this framework is not satisfactory for Markovian dynamics 
(5.19) on the Banach space £i(?£), although (5.19) is a strongly continuous 
^-automorphism even for unbounded generator L. 

To define a dynamics, which is dual to the Schrodinger picture (5.19), 
we recall that semigroup {£*}oo serves to calculate evolution of the normal 
states {w*(-)} t >o on observables A e Q3("H): 



u\A) = Tr n {C t {p)A) =: (£ 4 (p),A) , p € d(«) 



(5.22) 



Here 'B('H) denote the C*-algebra of bounded operators on the Hilbert space 
W and p is a density-matrix operator with the trace-norm ||p||e;i = 1. 
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Recall that by virtue of (5.22) the Banach space of bounded operators 
on H is topological^ dual of £i(H): 03(H) = (£i(H))*. This means that 
the map A 4 (-,A) is an isometric isomorphism of 03 (H) onto the set of 
linear continuous functionals ((£i(H))*, defined on the space £i(H). Then 
semi-norms generated by this duality 

{m p (A) ■■= KMMWiCW) > 03(H) , (5.23) 
define on the Banach space 03(H) the weak* -topology, which coincides with 
the operator c-weak topology, and for the operator norm of A one gets: 

II4|= sup (5.24) 



see e.g. |ReSiI| or [A JPIj . Duality (5.22) defines also the adjoint semigroup 
{£**}t>o on the dual space (£i(H))* = 03(H): 

(£ t (p),A) = (p,£ t *(A)) , pe€i(U) , Ae (£i(H))* . (5.25) 

In general, the adjoint semigroup is not strongly continuous on the dual 
Banach space (£i(H))*, although (5.25) and the strong continuity of (5.19) 
trivially imply the weak* -continuity of {£**}t>o on this space. 

Even thought {£**}t>o is not necessarily strongly continuous on 03(H), 
one can still associate with this semigroup a generator L in the weak*- 
topology: 

LA :=w* -lim^iC** A- A) , (5.26) 

with domain 

dom(Z) := {A G 03(H) : w* - hm Q -(C t *A - A) 3}. (5.27) 

It turns out that generator L is weak*-densely defined, closed operator, which 
coincides, LA = L*A, with the adjoint operator L* , i.e. 

dom(Z) = dom(L*) := (5.28) 

{A G 03(H) : 3 A' G 03(H) s.t. (p,A') = (L(p),A) for all p G dom(L)} . 

Although the weak*-topology is even weaker than the weak topology on the 
Banach space 03(H) the above arguments make legitime the characterization 
of semigroup {£' * :— e tL }t>a by a generator in a close similarity to strongly 
continuous case, see |BrRol| . Ch.3 for further details. The minimal price is 
that instead of the C*-algebra of bounded operators 03(H) one must consider 
this space endowed by a weaker, namely the weak* -topology. 

Recall that von Neumann algebra is a C*-algebra acting on H, con- 
taining identity operator and closed in the weak operator topology. It has 
enough room for the weak* -continuous semigroup {e tL }t>o- The Banach 
space 03(H) is example of a von Neumann algebra, while the Banach space 
<£i(H) is evidently not. 

Let SOT C 03(H) be a von Neumann algebra and the weak* -continuous 
for all A G 9JT map t H> r*(A) be a (semi)group of ^-automorphisms of 9Jt. 
Then the pair (9Jt, r*) is called a VK*-dynamical system. In our case r* = £**. 
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If the semigroup (5.19) is trace-preserving (5.9), then the adjoint semi- 
group (|5.25| is a unity-preserving (£**(!) = 11) contraction: ||£**(A)|| < \\A\\, 



see (5.24). This dual map inherits the property to be a completely positive 



semigroup, which was established for Markovian dynamics (5.19), see |Davl) 
and [Dav2 . By consequence, one gets for generator L* the analogue of the 
Kossakowski-Lindblad representation |FrVe[ lAlFaj : 



L* (A) =i[H,A] + l [A, V a ] + [V* ,A]V a ), A e <&{U) , 



(5.29) 



see (5.201, (5.21) and (5.28) 



Duality (5.25) is useful for control the state evolution £'(/?) and, in 
particular, for the proof of the t — ¥ oo limit p^ by calculation of this limit 
on observables. Since for any r > and p £ £i(H) , A £ Q5("H): 



Tr« (poo A) = lim Tr H (C T+t (p), A) = Tr n (C T ( Poo ),A) , 



(5.30) 



we conclude that p^ is £ r -invariant (steady) state. To elucidate topology of 
the density matrix convergence, recall that for the Kossakowski-Lindblad gen- 
erator (5.20), ( 5.21 1 dynamics of state is trace-preserving (5.9): ||£ t (p)||,>; 1 = 
1, and that (5.30) implies the weak- operator convergence C'(p) — > p^. Then 



(5.30) is equivalent to the trace-norm convergence of density matrices: 



lim \\C\p) - p.,, || 

t—too 







(5.31) 



see e.g. [Za|, Ch.2.4. 



A. 3 W* -Dynamics, Steady States and Return to Equilibrium. For open cavity 



(5.5) the choice of operators {VqIq in (5.29) is defined by (5.7). Since the 
cavity is a boson system, the natural (Fock) representation of the Canonical 
Commutation Relations (CCR) involves one-mode unbounded creation and 
annihilation operators b* and b acting on the Hubert space H = He- It 
is a boson Fock space H = 3b (C) over the one-dimensional subspace \) = 
{C 0}<cec (°f a Hilbert space) corresponding to this one photon mode 0. 

To avoid the problems with unbounded operators and to keep the evo- 
lution of observables in the space 03 (H) one considers the corresponding 
*-algebra of bounded Weyl operators I >rI!o2 in the form: 



2U(C) := <^ W(0 = exp 



V2 



((b + Cb*) 



C6C 



that verify the Weyl CCR-relations 

W(Ci)W(C 2 ) = e - lIm ^>/ 2 W(Ci + Ca) , 



(5.32) 



(5.33) 



Note that the family {W(()}^ & c is continuous on T-L in the strong operator 
sense, but it is not continuous in the C*-algebra topology: ||W(£) — 1|| = 2 
for any ( ^ 0. 
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By virtue of (5.7) and (5.21) one gets for the adjoint generator (5.29) 



of the one-mode open cavity with He = eb*b 
L* c JA)=i[eb*b,A] 



(5.34) 



-(<r-b*[A,b] + (T-[b* ,A]b + a+b[A,b*] + (T + [b,A]b*) , Ae<B{U) 



Then the adjoint semigroup equation: 

d t C t c %{A)=C t c %{L*cA A )) 



(5.35) 



} , m 



iet-{cr-— <r+)t/2 



allows to calculate evolution of the Weyl operators (5.32): A = W(Q, explic- 
itly: 

Cro= e " fl,(0 W (W)) • (5-36) 

Here 

4 <7_ — (7+ I 

(5.37) 

Since ||W(C) — 1|| = 2, the evolution (5.36) is not continuous in the 
C*-algebra topology, but it does i n the weak* -topology on the von Neumann 
algebra 2B(C) generated by (5.32) and the weak operator closure. Hence, the 
pair (2U(C),££* CT ) is W^-dynamical system, see [BrRoll lAlPl] . 

Note that by differentiating of C t c * a (W(C)) with respect to £ and £ one 
can calculate the evolution of polynomials of creation-annihilation operators 
in the weak*-topology. For example of the photon number operator N(t) := 
C t c * a (b*b), the formal evolution equation follows directly from (5.34) for A = 
b*b, and from ( |5.35| ): 

d t N(t) = -(cr_ - a+)N(t) + a+ , N(t = 0) = 6*6 

If a— > <J+ (leaking is stronger then pumping) , then one gets 



N(t) = e-^— a +^ 6*6+ —?± — (] 

<T_ — (7+ I 

By consequence, (5.39) formally implies 



(5.38) 
(5.39) 



lim C% (6*6) = 11 

t— >oo (J_ — (7 j 



and (5.36) and (5.371 yield 



lim C* c * a (W(0) = 1 exp{ 

t— >oo ' 4 



ICI 2 



<7_ 



} 



(5.40) 



(5.41) 



in the weak* -topology. 

Let the initial cavity state p be such that N(0) 



and (T_ > cr+ > 0. Then the limit (5.40) implies a nontrivial stationary 

(5.42) 



expectation value of the photon number: 



N(t) := lim Tr„ c (p£ £ *(6*6)) = lim Tr« c (£ £ (p) 6*6) 

t— too t— >oo 

= Tr «c (/°oo 6*6) = — — — , 
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in the limiting cavity state p x :— limt_ i . 00 Similarly, by (5.41 ) one gets: 



lim Tr« c (p£**(^(C)) = lim Tr« c (£ 4 (p) W(Q) 

t — Voo t— >oo 



(5.43) 



Tr nc ( Poo W(0)=eM- 



ICI' 1 



-} 



By (5.32) and (5.42), (5.43) we obtain that the limiting (steady) density 
matrix p^ of the open cavity corresponds to the one-mode boson equilibrium 
state: 

1 



Po 



-: p 



Pea 



(1 



-Pc* 



-ft av t b'b 



In — 

e a. 



(5.44) 



Therefore, if rr_ > a+ > and e > 0, the one-mode pumped leaky 
cavity evolves from any initial state verifying N(0) < oo, to the equilibrium 
Gibbs state with temperature 9 cav = l//3 cav entirely defined by the leaking- 
pumping intensities a^. For this limit of return to the thermal equilibrium 
one can distinguish two intuitively clear extreme cases. The zero-temperature 
case for zero pumping: a + = 0, and the infinite temperature case, when the 
pumping is not dominated by the leaking: a + j" cr_. In the first case the 



photon-number mean- value (5.42) is zero, whereas it is infinite in the second 
case. 

Note that the time and the pumping limits do not commute. Let the 
initial expectation of photons in the cavity N(0) < oo. Then applying the 
limit (7+ f <7_ to the expectation of ( |5.39| ) we obtain 



lim 



N(t) = N(0) + a + t 



(5.45) 



Therefore, one has a linear asymptotic increasing of the photon-number 
mean- value in the limit, when the leaking and pumping rates coincide: cr_ = 
o+ > 0. 

A.4 Completely Positive Quasi- Free Dynamics on the CCR Algebra. We 

need a more abstract description of dynamics £^,* (T on the CCR algebra 



2B(C) generated by (5.32). 



Notice that the complete positivity of the map C}^^ : 2U(C) -> 2U(C) 



follows from general properties of the Kossakowski-Lindblad generator (5.29). 



For the case of the open cavity this result follows directly from CCR-rclations 



(5.33) and the explicit result (5.36). Behind this result there are abstract 



observations due to |DVV11 lDVV2j and |EvLe2HVanj . 

Let 9) be a Hilbert space. Denote by T : fj —> S) a linear map and by 
: io — ► C a complex function. Consider the Weyl CCR(-f)) algebra generated 
by unitaries: 



W(Sj) := W(f) = exp 



V2 



W) + b*(f)) 



(5.46) 



fesi 



Here linear and anti-linear functions: / i— > b*(f) and / i— > b(f), are creation 
and annihilation operators in the boson Fock space $b ($)) over fj. Then the 
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Weyl CCR-relations take the form 

W(f)W(g) = e-* Im( ^^ 2 W(f + g) , f,g€S), (5.47) 

where (-,-)i5 i s the scalar product in Si. 

Recall now that a linear map T : 2B(f)) — > 2U(.f)) on the CCR(i^) algebra 
(5.461 is called quasi-free, if it has the form, |DVV1[ |LWV2| : 

T(W(f)) = *(/) W(T(f)) , feSj. (5.48) 

Then the unity-preserving quasi-free semigroup {T t } t > on the CCR(f)) al- 
gebra W(Sj) is defined in a similar way: 

T t (W(f)) := * t (/) W(r t (/)) , / G Si , (5.49) 

where * t=0 (/) - 1, r i=0 = 1 and * t (/ = 0) = 1, T t (f = 0) = . Note that 



the semigroup property of {T t } t >o and (5.491 imply 



T a+t (W(f)) = T s (T t (W(f))) 



(5.50) 



(5.51) 



= *.(r t (/))* t (/) w(r a (r t (/))) = * 5+t (/)w(r s+t (/)) . 

Then linear independence of the Weyl operators yields 

r s+t (f) = r a (r t (/)) and * 5+t (/) = *,(r t (/))* t (/) . 

Hence, {r t } t > is in turn a semigroup on fj. 

Let t i-> **(/) be continuous and for each / S f) be differentiable at 
t = +0: 

$i(/):= t lim9 ( $ ( (/), (5.52) 



such that the function 1 i— > <3>Q(r t (/)) be bounded. Since semigroup (5.491 is 
unity-preserving, then (5.51 ) and W(f = 0) = 11 imply that for any / £ S) 



* t (/)=exp / dr%(T r (f)) 



(5.53) 



where %{0) = and %(-f) = %{f), see |DVV2j 

Note that for the particular case: Sj = C and for identification of (|5.32|) 



with (5.46), we reproduce the results (5.361 and (5.371 in A. 3 for T t =C? C a . 
To this end one has to use (5.51 1, ( 5.52[ ) and to put 



IMC) :=(e i£t -( ff - ff +"/ 2 and *' (C) 



|C| 2 



(o-_+cr+). (5.54) 



Then (5.53) gives * t (C) 



-n*(0 



where flt(C) is defined by (5.37) 



A. 5 Quasi- Free States on the CCR Algebra. Since the linear hull of the the 

Weyl operators is dense in CCR(f)), any state lj(-) is uniquely determined by 
its values taken on {W(/)}/efl- Therefore a state uj is completely defined by 
its characteristic functional 

fy 3 f h+ u>{W{f)), 

A state to is called regular if the function a > ui(W(af)) is continuous for all 
/ € Sj. Characteristic functionals of regular states on CCR(i}) are character- 
ized by Araki and Segal [AJPI, BrRol] in the following theorem. 
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Theorem 5.2. A map $) 9 / >— > o>(W^/)) € C is £/ie characteristic functional 
of a regular state u) on CCR(S) ) if and only if 

1. lo{W{Q)) = 1. 

2. The function a <— > uj(W(af)) is continuous for all f 6 $). 

3. For any integer n > 2, aZZ /i, G .fj anrf aZZ z%, ...z n G C one has 

n 

J2 u{W{f 3 - f k ))e- llm ^<^/ 2 z]z k > 0. 

j,k=l 

We remind that the state ui r . s (-) is called quasi-free | Ver] . if 

Wr,.(W r (/))=exp{ir(/)-i «(/,/)} , /Gij. (5.55) 

Here r is a linear functional on S), whereas s is a non- negative (closable) 
sesquilinear form on fj x Sj, that verifies 

- A |Im(/, | 2 < s(/, /) s( ff , 3 ) , f,g€f>, (5.56) 

to ensure the positivity of this state (Ver]. By the Araki-Segal theorem the 
quasi-free states are regular and analytic, verifying the equations: 

r(/)=w r ,.(*(/)) and s(f, /) = ^ s ($(/) 2 ) - co,, s (Hf)f , (5.57) 

where $(/) := (b(f) + b*(f))/V2, see, e.g. [VeR IBrRo2] . 

If the state (5.55) is gauge-invariant: uj ryS (W(f)) = u) riS (W(e l v /)), 
then r(-) = and we denote this state by lu s (-) := ui r —Q tS (-). By virtue of 
( 5.43[ ) the limiting (steady) state (5.44) of the open cavity is gauge invariant 
and quasi-free with: 

r(0 = o , s(c C) = 4r = Tr « c (/ ca " $ (0 2 ) - C e £ = C , (5.58) 

l er_ — cr_)_ 

Moreover it is the equilibrium, i.e. (/3 ca „)-KMS state with respect to the 
Hamiltonian dynamics r* generated by He- 

It is worth to note that a priori there is no evidence that the steady 
state of the open cavity (5.5 1 must be either quasi-free, or the Gibbs equilib- 
rium state. For discussion of the theory of non- equilibrium quasi-free steady 
states we suggest a very complete account in |AJPI] - |AJPIIIj . There one 
can find curious examples of non-equilibrium (non-KMS) quasi-free states, 
which allow a certain informal Gibbs description via long-range many-body 
interactions [AschPi] . 

Note that by definitions (5.49) and (5.55) the quasi-free dynamics maps 
the quasi-free states into the states: 

uj r , s (T t (W(f))) = Mf) Wr,.(W(T t (/))) - *,(/) oj ruSt (W(f)) , (5.59) 

where r t (f) := r(T t (f)) and s t (f,f) := s(T t (f),T t (f)). In general the states 
(5.59) are not quasi- free. 



48 B. Nachtergaele, A. Vershynina and V. A. Zagrebnov 



Let us consider the case of Gaussian quasi-free dynamics [Vanj . when 
the semigroup {r^oo defined on $j by 



T t :=exp{itH- it(E_ - £+)} 



(5.60) 



here H is a self-adjoint operator and Sip are bounded positive operators on 
Sj such that S_ > S + > 0, and (5.52) is a bilinear form 



*o(/) '=-\(f,Rf) , /efl, 



(5.61) 



defined by a positive bounded operator R > S_. Then by virtue of (5.53), 
( 5.55 1 and ( 5.59 ) dynamics T t maps initial quasi-free state w riS into quasi-free 
state u) rtt j t with 



*t(f,f) :=*(Tt(/),r t (/)) + 



dr(r T (/),i?r T (/)) 



(5.62) 



In the particular case of Hamiltonian dynamics (Szp = and R = 0) the 
quasi-free map (5.49) is the group of Bogoliubov automorphisms on W($j): 

T t (W(f)) = W(e ltH f) , feSj. 



(5.63) 

Automorphism (5.63) is the simplest quasi-free dynamics, which corresponds 
to the unitary one-particle evolution generated by the Hamiltonian H. 

For the example of the open cavity, when Sj = 
establish that in this case one has to put H 
( 5~60| and ( |5~61~] ). Then ( pTo^ ) yields 

^(C,0 = S (r t (C),r t (C)) + i|CI 2 a -^ ±1± (i- e -(—^*) . (5.64) 

z cr_ — tr_i_ 



Sir = <Jzr, R 



we use (5.54) to 
cr_ + cr + in 



This means that dynamics generated by (5.34) is quasi-free: C\^ a {W(()) 
T t (W(C)), see (5.54), and that it preserves the quasi- free states. 
Since lim^oo r t (£) = 0, by (5.64) we obtain that 

Id 2 a-^~ 



lim Wr , s (£^(W(C))) = exp{-J 



(5.65) 



for any initial quasi-free state io r 



In fact we established in (5.43) that for t — > oo the quasi-free dynamics 
C t c * a transforms any initial state p with finite expectation of photon number 
into the limit state, which is quasi-free and Gibbs, see (5.43), (5.44) and 
(5765 1. 



Non-equilibrium states of a photon cavity 



49 



References 

[Al] R. Alicki, Invitation to Quantum Dynamical Semigroups, pp 239-264, in 
Dynamics of Dissipation P. Garbaczewski and R. Olkiewicz (Eds.), Lec- 
ture Notes in Physics 597, Springer Verlag, Berlin-Heidelberg 2002. 

[AlFa] R. Alicki and M. Fannes, Quantum Dynamical Systems, Oxford University 
Prass, Oxford 2001. 

[AschPi] W. Aschbacher and C.-A. Pillet, N on- Equilibrium Steady States of the XY 
Cham, J. Stat. Phys. 112 (2003) 1153-1175. 

[APa] S. Attal and Y. Pautrat, From repeated to continuous quantum interac- 
tions, Ann. Henri Poincare 7 (2006) 59-104. 

[AJ1] S. Attal and A. Joye, Weak coupling and continuous limits for repeated 
quantum interactions, J. Stat. Phys. 126 (2007) 1241-1283. 

[AJ2] S. Attal and A. Joye, The Langevin equation for quantum heat bath, J. 
Func. Anal. 247 (2007) 253-288. 

[AJPI] Open Quantum Systems I, The Hamiltonian Approach, S. Attal, A. Joye, 
C.-A. Pillet (Eds.), Lecture Notes in Mathematics 1880, Springer- Verlag, 
Berlin-Heidelberg 2006. 

[AJPII] Open Quantum Systems II, The Markovian Approach, S. Attal, A. Joye, 
C.-A. Pillet (Eds.), Lecture Notes in Mathematics 1881, Springer- Verlag, 
Berlin-Heidelberg 2006. 

[AJPIII] Open Quantum Systems III, Recent Developements, S. Attal, A. Joye, 
C.-A. Pillet (Eds.), Lecture Notes in Mathematics 1882, Springer- Verlag, 
Berlin-Heidelberg 2006. 

[BrRol] O. Bratteli and D.W. Robinson, Operator Algebras and Quantum Statis- 
tical Mechanics,vo\.l, Springer- Verlag, Berlin 1979. 

[BrRo2] O. Bratteli and D.W. Robinson, Operator Algebras and Quantum Statis- 
tical Mechanics, vol. 2, Springer- Verlag (2nd Edt), Berlin 1997. 

[BJM1] L. Bruneau, A. Joye, and M. Merkli, Asymptotics of repeated interaction 
quantum systems, J. Func. Anal. 239 (2006) 310-344. 

[BJM2] L. Bruneau, A. Joye, and M. Merkli, Random repeated interaction quantum 
systems, Commun. Math. Phys. 284 (2008) 553-581. 

[BJM3] L. Bruneau, A. Joye, and M. Merkli, Repeated and continuous interactions 
in open quantum systems, Ann. Henri Poincare 10 (2010) 1251-1284. 

[BPi] L. Bruneau and C.-A. Pillet, Thermal relaxation of a QED cavity, J. Stat. 
Phys. 134 (2009) 1071-1095. 

[ChFa] A.M. Chebotarev and F. Fagnola, Sufficient conditions for conservativity 
of minimal quantum dynamical semigroups, J. Func. Anal. 153 (1998) 382- 
404. 

[Davl] E.B. Davies, Quantum dynamical semigroups and the neutron diffusion 

equation, Rep. Math. Phys. 11 (1977) 169-188. 
[Dav2] E.B. Davies, Generators of dynamical semigroups, J. Func. Anal. 34 (1979) 

421-432. 

[DVV1] B. Demoen, P. Vanheuverzwijn, A. Verbeure, Completely Positive Maps 
on the CCR-Algebra, Lett. Math. Phys. 2 (1977) 161-166. 

[DVV2] B. Demoen, P. Vanheuverzwijn, A. Verbeure, Completely Positive Quasi- 
Free Maps on the CCR-Algebra, Rep. Math. Phys. 15 (1979) 27-39. 



50 



B. Nachtergaele, A. Vershynina and V. A. Zagrebnov 



[EvLel] D.E. Evans and J.T. Lewis, Dilations of irreversible evolutions in algebraic 
quantum theory, Commun. Dublin Inst. Adv. Studies, Ser. A32, 1977. 

[EvLe2] D.E. Evans and J.T. Lewis, Some semigroups of completely positive maps 
on the CCR algebra, J. Funct. Anal. 26 (1977) 369-377. 

[FNW] M. Fannes, B. Nachtergaele, and R.F. Werner, Comm. Math. Phys 144 
(1992) 443-490. 

[FJMa] P. Filipowicz, J. Javanainen, and P. Meystre, Quantum and semiclassical 
steady states of a kicked cavity mode. J. Opt. Soc. Am. B 3 (1986) 906-910. 

[FJMb] P. Filipowicz, J. Javanainen, and P. Meystre, Theory of a microscopic 
maser, Phys. Rev. A 34 (1986) 3077-3087. 

[FaRe] F. Fagnola and R. Rebolledo, Entropy Production for Quantum Qynamical 
Semigroups, arXiv:1212.1366vl [math-phys] 6 Dec 2012. 

[FrVe] A. Frigerio and M. Verri, Long-time asymptotic properties of dynamical 
semigroups on W* -algebras, Math. Z. 180 (1982) 275-286. 

[Hoi] A.S. Holevo, On dissapative stochastic equations in a Hilbert space, Prob. 
Theor. Rel. Fields 104 (1996) 483-500. 

[MWM] D. Meschede, H. Walther, and G. Miiller, One-atom maser, Phys. Rev. 
Lett. 54 (1985) 551-554. 

[NVZ] B. Nachtergaele, A. Vershynina, and V.A. Zagrebnov, Lieb-Robinson 
Bounds and Existence of the Thermodynamic Limit for a Class of Ir- 
reversible Quantum Dynamics, Contemporary Mathematics (AMS) 552 
(2011) 161-175. 

[NZ] H. Neidhardt and V.A. Zagrebnov, Linear non- autonomous Cauchy prob- 
lems and evolution semigroups, Adv. Differential Equations 14 (2009), 289- 
340. 

[ReSil] M. Reed and B. Simon, Methods of Modern Mathematical Physics I. Func- 
tional Analysis, Academic Press, New York 1972. 

[Si] K.B. Sinha, Quantum Dynamical Semigroups, pp. 161-169, in Operator 
Theory: Advances and Applications 70, Birkhauser-Verlag, Basel 1994. 

[Van] P. Vanheuverzwijn, Generators for quasi-free completely positive semi- 
groups, Ann. Inst. Henri Poincare XXIX (1978) 123-138. 

[Ver] A. Verbeure, Many-Body Boson Systems, Half a Century Later, Springer- 
Verlag, Berlin- Heidelberg 2011. 

[VWZ] P.-A. Vuillermot, W.F. Wreszinski, and V.A. Zagrebnov, A general 
Trotter-Kato formula for a class of evolution operators. J. Funct. Anal. 
257 (2009) 2246-2290. 

[Za] V.A. Zagrebnov, Topics in the Theory of Gibbs Semigroups, Leuven 
Notes in Mathematical and Theoretical Physics - Series A: Mathemati- 
cal Physics, Volume 10, Leuven University Press 2003. 

Bruno Nachtergaele 
Department of Mathematics 
University of California, Davis 
Davis, CA 95616 
USA 

e-mail: bxn@math . ucdavis . edu 



Non-equilibrium states of a photon cavity 



Anna Vcrshynina 
Department of Mathematics 
University of California, Davis 
Davis, CA 95616 
USA 

e-mail: annavershyninaOgmail . com 

Valentin A. Zagrebnov 

Laboratoire d" Analyse, Topologie, Probabilites - UMR 7353 
CMI - Technopole Chateau- Gombert 

39, rue F. Joliot Curie, 13453 Marseille Cedex 13, France 
and 

Departement de Mathematiques, Aix-Marseille Universite 
Luminy - Case 901, Marseille 13288, Cedex 09, France 
e-mail: Valentin. Zagrebnov@latp .univ-mrs . f r 
Valentin. Zagrebnov@univ-amu. f r 



